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Although  turbulent  boundary  layers  are  very  common  flows  in  engineering,  the 
physical  mechanisms  responsible  for  the  production  and  sustenance  of  turbulence  are 
still  ill-defined.  Largely  due  to  visual  studies,  the  traditional  view  of  turbulence  as  a 
random  field  has  evolved  in  the  past  twenty  years  into  a  more  mechanistic  description 
of  the  turbulent  processes  through  the  identification  of  so-called  coherent  structures. 
Among  the  various  structures  identified  in  the  context  of  turbulent  boundary  layers  are 
the  low-speed  streaks  (Schraub  and  Kline,  1965;  Kline  et  al.,  1967;  Blackwelder,  1978 
and  others),  the  hairpin  vortices  (Kim  et  al.,  1971;  Smith,  1978;  Head  and  Bandyopha- 
dyay,  1981,  Moin  and  Kim,  1985  and  others)  and  the  large-scale  outer  structures 
(Blackwelder  and  Kovaznay,  1972;  Brown  and  Thomas,  1977;  Chen  and  Blackwelder, 
1978  and  others).  It  has  also  been  recognized  that  most  of  the  turbulence  production 
occurs  near  the  wall  (Klebanoff,  1954)  during  very  localized  random  events  (in  space 
and  time)  called  bursts  (Kline  et  al,  1967;  Kim  et  al,  1971;  Blackwelder  and  Kaplan, 
1976  and  others).  Due  to  the  complexity  of  these  unsteady  three-dimensional 
processes,  and  the  lack  of  universally  accepted  definition  of  their  spatio-temporal  struc¬ 
ture,  most  of  the  evidence  gathered  from  measurements  or  visual  studies  has  been 
incomplete  and  often  hard  to  meld  into  a  unified  view  of  the  boundary  layer. 

Recent  advances  in  computational  and  instrumentation  means  have  triggered 
new  approaches  to  improve  our  understanding  of  boundary  layers:  first,  the  possibility 
of  undertaking  full  Navier-Stokes  simulations  for  simple  turbulent  flows  (Moin  and 
Kim,  1982,  etc...);  second,  the  possibility  of  processing  efficiently  large  images  of  flow 
fields  to  obtain  quantitative  information  about  their  instantaneous  structures  (Brodkey, 
1986  and  others);  and  third,  the  possibility  of  acquiring  large  amounts  of  experimental 
data  to  map  (at  least  in  the  statistical  sense)  the  three-dimensional  flow  field  (Guezen- 
nec,  1985  and  others).  It  is  only  through  the  convergence  of  these  efforts  that  the 
dynamics  of  the  boundary  layers  will  be  fully  elucidated,  and  that  previous  findings 
will  be  integrated  into  a  more  mechanistic  description. 

It  is  very  important  from  a  technological  standpoint  to  understand  the  physical 
processes  in  turbulent  boundary  layers,  especially  to  achieve  desired  characteristics  by 
modifying  their  structure.  For  example,  the  recent  results  in  drag  reduction  through 
suppression  of  large  scales  (Corke  et  al,  1982;  Plesniak  et  al.,  1985  and  others)  could 
be  further  optimized  if  the  mechanisms  at  play  were  better  understood.  Other  potential 
applications  include  heat  and  mass  transfer  modification  and  boundary  layer  noise 
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reduction.  The  ubiquity  of  these  wall-bounded  turbulent  flows  in  technological  appli¬ 
cations  justifies  the  need  for  a  more  fundamental  understanding  of  the  underlying  phy¬ 
sics. 


An  other  important  aspect  is  to  bridge  the  gap  between  transitional  and  turbulent 
boundary  layers.  Recent  studies  seem  to  indicate  that  similar  structures  and  physics 
may  be  present  in  both  flows  (Kim  and  Moin,  1986;  Williams,  1985  and  others),  but 
no  direct  attempt  at  tracing  these  structures  over  a  range  of  Reynolds  number  in  both 
regimes  has  yet  been  undertaken.  It  is  clear  that  if  the  same  physical  processes  are 
responsible  for  the  breakdown  of  transitional  boundary  layers  and  the  turbulence  pro¬ 
duction  in  turbulent  boundary  layers,  major  advances  in  the  modeling  of  such  flows 
way  can  be  achieved  by  a  unified  approach. 

Nagib  and  Guezennec  (1986)  performed  the  first  attempt  to  map  experimentally 
the  three-dimensional  structure  of  some  turbulence  producing  events.  Their  results 
indicate  the  presence  of  two  counter-rotating  roller  structures  with  a  scale  of  the  order 
of  half  the  boundary  layer  thickness.  Based  on  these  results,  a  conjectured  model  of 
the  turbulent  boundary  layer  was  proposed.  However,  a  number  of  questions  remained 
unanswered.  The  origin  of  those  rollers  was  not  clearly  established,  nor  was  their  evo¬ 
lution.  In  addition,  the  mapped  events  were  only  reconstructed  for  two  classes  of 
events  due  to  the  particular  detection  technique  used  (wall-shear).  Although  they  con¬ 
clusively  proved  that  the  events  detected  by  this  technique  were  second  and  fourth 
quadrant  Reynolds  stress  events,  a  broader  class  of  events  than  the  conventional  qua¬ 
drant  technique  may  have  been  averaged,  which  could  have  led  to  some  smearing  of 
the  spatial  maps.  This  was  later  shown  to  be  true  particularly  for  the  second  quadrant 
motions  by  Guezennec,  Piomelli  and  Kim  (1987).  Moreover,  due  to  e  xperimental 
limitations,  the  transverse  component  of  velocity  was  not  measured,  but  inferred  from 
continuity  assuming  symmetry  about  the  center  plane.  This  symmetry  which  reflects 
the  homogeneity  of  the  statistics  in  the  spanwise  direction  does  not  need  to  be  true 
instantaneously,  but  only  on  the  average.  In  addition,  it  has  been  speculated  that  that 
instantaneous  stmctures  may  come  isolated  rather  than  in  pairs,  as  implied  by  statisti¬ 
cal  averages.  Guezennec,  Piomelli  and  Kim  (1987)  and  Choi  and  Guezennec  (1990) 
quantified  this  asymmetry  and  found  that  instantaneous  structures  are  more  likely  to  be 
asymmetric.  Based  on  this,  they  designed  new  ensemble  averaging  procedures  to 
respect  this  marked  asymmetry. 

More  recently,  Guezennec  (1988)  reprocessed  his  database  using  the  stochastic 
estimation  technique  (Adrian,  1979;  Adrian  and  Tung,  1980,...).  The  results  obtained 
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by  this  technique  compared  very  well  with  those  obtained  by  conventional  ensemble 
averaging.  The  stochastic  estimation  offers  many  advantages  over  conventional  pro¬ 
cessing.  In  particular,  it  allows  to  examine  any  number  of  "conditions"  selected  a  pos¬ 
teriori  without  reprocessing  die  raw  data  base,  since  the  procedure  relies  upon  uncon¬ 
ditional  statistics  (two-point  spatial  correlations).  This  represent  a  real  computational 
advantage  when  dealing  with  large,  three-dimensional  data  bases.  Time  evolution  for  a 
given  condition  can  be  obtained  by  a  simple  extension  of  the  technique  and  the  use  of 
space-time  correlation  functions  instead  of  spatial  correlations.  Due  to  the  uncondi¬ 
tional  nature  of  the  underlying  statistics,  the  estimates  of  the  conditional  averages  also 
tend  to  be  smoother  than  those  obtained  by  conventional  means  (with  limited  number 
of  realizations). 

It  has  been  observed  that  single  point  conditional  averages  (conventional  or  else) 
tend  to  yield  structures  which  scale  with  the  integral  scale  of  the  flow.  One  way 
around  this  problem  is  to  perform  multi-point  conditional  averages,  i.e.  specifying  con¬ 
ditions  at  more  than  one  point.  This  allows  to  examine  events  with  a  specific  length 
scale  governed  by  the  separation  between  the  two  or  more  points  at  which  conditions 
are  specified.  However,  such  conditional  measurements  are  hard  to  obtain  in  practice 
by  ensemble  averaging.  On  the  other  hand,  the  stochastic  estimation  is  easily  general¬ 
ized  to  estimate  multi-point  conditional  averages.  Preliminary  applications  of  two-point 
conditional  averages  have  been  performed  by  Adrian  and  Moin  (1987)  and  Guezennec 
(1988)  to  examine  events  characterized  by  a  sharp  interface  such  as  VITA  events 
corresponding  to  the  rapid  interaction  between  second  and  fourth  quadrant  motions. 
The  multi-point  conditional  averages  is  a  very  valuable  new  tool  to  study  "realistic" 
average  structures  which  bear  some  resemblance  to  the  instantaneous  ones  (see  Adrian, 
Moser  and  Moin,  1987). 

II.  DESCRIPTION  OF  EXPERIMENTS 

The  experimental  research  performed  was  aimed  at  gaining  a  better  understand¬ 
ing  of  the  structure  of  turbulent  boundary  layers  and  the  mechanisms  responsible  for 
the  sustenance  of  turbulence.  In  addition,  the  relationship  between  structures  observed 
in  turbulent  and  transitioning  boundary  layers  were  to  be  examined.  These  objectives 
were  addressed  by  performing  three-dimensional  mappings  of  the  flow  field  at  various 
Reynolds  numbers  ranging  from  the  transition  regime  to  the  fully  turbulent  one. 

One  of  the  key  aspects  of  this  study  was  not  only  to  extract  the  important  flow 
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structures,  but  also  to  study  their  evolution.  Due  to  the  ease  of  estimating  the  condi¬ 
tional  flow  fields  for  a  large  number  of  conditions  by  means  of  the  stochastic  estima¬ 
tion,  the  evolution  of  the  structures  can  then  be  obtained  by  reconstructing  a  sequence 
of  flow  fields  corresponding  to  a  real  (measured)  time  sequence  of  conditions  (pseudo- 
dynamics).  Essentially,  this  procedure  allows  to  examine  the  average  evolution  of 
structures  corresponding  to  real  sequences  of  events.  This  approach  should  enchance 
our  fundamental  understanding  of  the  dynamics  of  the  turbulence,  while  eliminating 
some  of  the  unnecessary  details  which  contributed  to  the  event-to-event  variability. 
The  realism  of  such  pseudo-dynamic  reconstruction  is  particularly  enhanced  by  using 
multi-point  conditional  estimates  as  described  earlier. 

Based  on  the  limitations  of  the  previous  work,  it  is  felt  that  all  three  components 
of  velocities  had  to  be  measured  both  at  the  detection  probe  and  at  the  moving  probe. 
The  penalty  paid  for  having  a  intrusive  detection  probe  (as  opposed  to  flush-mounted 
shear  sensors)  is  outweighed  by  the  benefits  of  the  full  determination  of  the  velocity 
vector  and  allows  to  discriminate  a  broader  range  of  motions,  in  particular  the  more 
likely  asymmetric  structures.  The  experiments  consisted  of  measuring  the  full  three- 
dimensional  two-point  correlation  tensors  for  all  velocity  components  in  the  boundary 
layer  at  various  Reynolds  numbers  and  various  reference  heights.  The  measurements 
were  performed  using  a  three-component  miniature  hot-wire  probe  located  at  a  fixed 
position  above  the  wall  and  acquiring  data  simultaneously  with  another  three- 
component  probe  located  at  various  positions  around  the  fixed  probe.  The  movable 
probe  scanned  a  full  three-dimensional  grid  around  the  fixed  probe.  All  data  was 
acquired  digitally  and  calibrated  on-line.  Due  to  the  large  amount  of  data  required  to 
obtain  converging  statistics,  those  statistics  (two-point  spatial  correlations  of  second 
and  third  order)  were  also  computed  on-line  and  stored  for  each  grid  point.  Nonenthe- 
less,  the  raw  hot-wire  signals  were  archived  on  magnetic  tapes  after  calibration  to 
allow  future  reprocessing  of  the  data  by  other  means. 

Only  the  linear  estimation  technique  mentioned  earlier  was  used  to  process  the 
correlations  and  obtain  the  spatial  arrays  of  coefficients  required  to  estimate  the  aver¬ 
age  flow  field  for  any  given  condition.  While  the  quadratic  estimates  which  involve 
third  order  statistics  have  been  found  somewhat  better  to  examine  flow  fields  character¬ 
ized  by  a  high  skewness  (Guezennec,  1988),  this  advantage  disappear  when  multipoint 
linear  stochastic  estimates  are  used. 

As  mentioned  earlier,  the  evolution  of  the  important  structures  was  investigated 
by  measuring  real  time  sequences  of  conditions  at  the  fixed  probe  and  using  those  to 
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generate  dynamic  sequences  of  the  corresponding  averaged  flow  fields.  This  novel 
procedure  which  we  call  pseudo-dynamics  allows  to  obtain  a  animated  sequence  of 
kinematic  reconstruction  of  the  flow  fields  consistent  with  a  dynamic  (measured)  time 
sequence  of  conditions.  Animated  display  and  movies  of  such  sequences  of  events  are 
currently  being  generated  to  facilitate  the  interpretation  of  such  data. 

Since  the  mappings  of  the  boundary  layer  were  performed  for  all  velocity 
components,  it  is  more  efficient  to  measure  simultaneously  all  three  components  of  the 
velocity.  However,  the  large  amount  of  data  required  for  the  three-dimensional  map¬ 
pings  necessitate  veiy  efficient  calibration  and  velocity  extraction  procedures  for  the 
triple  hot-wire  probes.  Such  probes  are  known  to  be  very  difficult  to  calibrate  and  use. 
The  extraction  of  the  velocity  components  from  the  three  anemometer  voltages  is  a 
complex  task  involving  solving  a  stiff  system  of  non-linear  equations,  for  every  velo¬ 
city  triplet.  Hence,  a  considerable  amount  of  time  was  devoted  at  the  beginning  of  the 
project  to  establish  efficient  calibration  procedures  and  develop  a  novel  table  look-up 
technique  for  the  velocity  extraction.  The  details  of  such  techniques  can  be  found  in 
T.  J.  Gieseke’s  Masters  Thesis  (1990)  and  are  summarized  in  a  paper  presented  at  the 
Symposium  on  Turbulence  at  the  University  of  Rolla,  Misourri  in  September  1990  and 
accepted  for  publication  in  Experiments  in  Fluids.  A  copy  of  this  paper  is  included  in 
Appendix  A. 

Two  case  studies  were  conducted.  One  corresponded  to  the  experiments 
described  originally  in  the  proposal,  namely  a  mapping  of  a  regular  boundary  layer 
with  respect  to  a  probe  located  in  the  near  wall  region  and  performed  at  various  Rey¬ 
nolds  numbers.  The  other  one  was  a  comparative  study  of  the  large  scale  structures  in 
a  regular  turbulent  boundary  layer  and  the  changes  imposed  on  such  structures  in  the 
near  field  of  flow  manipulators,  namely  Large  Eddy  Break  Up  Devices. 

Case  1:  Wall  Laver  Structure  as  a  Function  of  Reynolds  Number 

For  this  experiment,  simultaneous  readings  of  velocity  were  taken  with 
two  triple  hot-wire  probes  placed  in  the  turbulent  boundary  layer.  One  was  fixed  in 
space  being  attached  to  the  measurement  plug  (the  reference  probe)  and  one  was 
attached  to  the  traversing  mechanism  (the  grid  or  mapping  probe).  The  mapping  probe 
was  moved  over  a  grid  of  locations.  These  experiments  were  performed  for  Reynolds 
numbers  of  3700,  2540,  and  960.  For  each  R  q,  the  experiment  was  repeated  for  two 
reference  probe  locations,  each  location  depending  on  Reynolds  number  and  experi¬ 
mental  grid  specifications.  A  triple  wire  probe  attached  to  the  traversing  mechanism,  a 
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triple  wire  probe  was  mounted  beneath  the  test  plate  with  the  hot-wire  array  and  body 
protruding  through  a  small  hole  in  the  measurement  plug  as  shown  in  Figure  1.  This 
probe  could  be  moved  vertically,  toward  or  away  from  the  test  plate,  by  sliding  it  in  a 
specially  designed  mount.  This  allowed  the  probe  to  be  moved  out  of  the  boundary 
layer  for  calibration  or  to  any  vertical  location  for  use  as  a  reference  probe  in  the  con¬ 
ditional  eddy  estimation.  Only  one  reference  probe  was  mounted  in  the  test  section  at 
a  given  time.  Additional  reference  points  were  added  by  repeating  the  experiments 
and  combining  the  results  in  processing. 

The  probe  mounted  to  the  traversing  mechanism  was  type  AHWP-3-2- 
OSU  (Probe  1).  The  stationary  probe  was  type  AHWP-3-1-OSU  (Probe  2).  Physical 
limitations  of  the  arrangement  came  from  the  size  of  the  probes.  The  delicate  structure 
of  the  hot-wire  arrays  required  that  extreme  care  be  taken  not  to  allow  the  probes  to 
come  into  contact  with  each  other  which  would  cause  damage  to  the  probes.  A  tele¬ 
scope  mounted  on  a  manual  positioning  mechanism  with  a  range  of  0  to  1.3  cm  in 
each  direction  was  used  to  accurately  measure  the  dimensions  of  each  probe  (Figure 
2).  This  allowed  definition  of  a  safe  range  for  the  moving  mapping  probe  around  the 
stationary  probe.  The  motion  of  the  probes  was  most  limited  in  the  z  direction  where 
the  probe  bodies  interfered  with  each  other.  A  computer  controlled  the  probe  motion 
using  an  algorithm  designed  to  avoid  objects.  These  objects  were  defined  in  space  by 
the  user  and  the  control  of  probe  motion  prevented  probe  contact  while  allowing 
minimum  probe  separation. 

The  grid  over  which  the  moving  probe  was  traversed  varied  with  test  Rey¬ 
nolds  number.  For  all  Reynolds  number  cases,  the  reference  probe  was  positioned  as 
close  to  the  wall  as  possible,  0.81  mm  when  measured  from  the  wall  to  the  center  of 
the  probe  measurement  volume.  This  brings  the  lowest  probe  prongs  nearly  in  contact 
with  the  wall  at  a  height  of  0.43  mm.  The  grid  for  each  case  is  described  by  the  infor¬ 
mation:  grid  comer  location  {x^^,yh,Zf^\  grid  size  {Sx,Sy,s^),  number  of  grid  points 
{n^,ny,n^),  reference  probe  location  {Xp,yp,Zp),  and  bounds  of  the  box  around  the 
reference  probe  to  be  avoided  {x^^-Xg<x<Xh+Xg,  0<y<yg,  Zf^-z^<z<z^^+Zg)  as  shown  in 
Figure  3.  Grid  dimensions  were  selected  to  capture  the  behavior  of  large  scale  eddies 
in  space  and  time  as  found  in  the  experiments  by  Choi  and  Guezennec  (1989).  Digital 
sampling  of  the  anemometer  readings  was  conducted  with  a  At  of  two  viscous  time 
units.  This  frequency  varies  roughly  like  the  square  of  the  Reynolds  number  as  the 
turbulent  time  scales  do.  The  number  of  samples  chosen  for  statistical  convergence  had 
to  be  reduced  from  20,000  to  10,000  for  the  lowest  Reynolds  number  case  in  order  to 
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Figure  1,  Schematic  Diagram  of  Probe  Mounting 
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Table  1.  Mapping  Grid  Characteristics 
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allow  the  experiment  to  be  completed  in  a  reasonable  time.  These  variables  for  all 
experiments  are  shown  in  Table  1  with  values  in  wall  units.  Owing  to  the  usual  sym¬ 
metric  structure  of  the  correlation  tensor  in  the  spanwise  direction,  the  measurements 
were  only  conducted  on  one  side  of  the  reference  probe. 

With  the  grid  definition  based  on  the  wall  variables  and  the  physical  limita¬ 
tions,  the  experiments  proceeded  in  a  straightforward  fashion.  Both  the  reference  and 
the  moving  probes  were  first  calibrated  to  determine  their  sensitivities  and  the  informa¬ 
tion  was  stored  as  a  set  of  coefficients  from  the  polynomial  that  best  fit  measured  vol¬ 
tages  and  freestream  velocities.  The  reference  probe  was  then  moved  into  position  in 
the  boundary  layer  and  the  mapping  probe  moved  to  a  known  location  near  the  refer¬ 
ence  probe.  With  the  free-stream  velocity  set  to  the  desired  value  computer  was 
given  control  of  the  experiment.  The  mapping  probe  traveled  from  point  to  point  on 
the  grid  where  the  record  of  discrete  voltages  output  from  each  of  the  six  anemometers 
was  stored  on  disk.  No  processing  of  the  readings  was  conducted  during  the  experi¬ 
ments.  The  raw  data  and  calibration  information  was  transfered  to  an  alternate  com¬ 
puter  system  where  calculations  could  be  conducted  more  rapidly. 

All  the  processing  took  place  on  a  DEC-3100  workstation  rather  than  on  the 
Masscomp  system  used  to  acquire  the  data.  This  was  primarily  done  because  the 
Masscomp,  although  ideal  for  data  acquisition,  is  much  slower  computationally  than 
the  DEC.  The  DEC,  however,  is  not  equiped  with  a  large  data  storage  system  like  the 
Masscomp  350  MB  hard  disk  drive  which  is  capable  of  storing  data  from  several 
experiments  before  a  tape  backup  is  required.  Because  of  the  limitations  of  both  sys¬ 
tems,  they  had  to  be  used  in  tandem.  Data  was  stored  on  the  large  disk  and  transfered 
a  parcel  at  a  time  to  the  DEC  for  processing.  The  data  parcels  were  organized  accord¬ 
ing  to  experimental  grid  location  and  contain  all  data  acquired  while  the  mapping 
probe  was  at  that  location.  Each  parcel  of  data  consisted  of  six  sampled  traces,  each 
trace  containing  10,000  or  20,000  data  points. 

Using  the  sensitivity  calibration  and  angular  calibration  tables  generated  for  the 
experiment,  the  voltage  readings  were  transformed  into  discrete  velocity  vector  traces 
at  the  reference  and  grid  locations.  After  the  mean  velocity  components  were 
removed,  correlations  between  the  fluctuating  velocity  components  were  calculated  to 
generate  the  stochastic  estimation  matrices  at  each  location.  In  addition  to  calculating 
the  correlations  between  simultaneously  occurring  velocities,  the  correlations  were  also 
computed  between  velocities  at  the  reference  position  and  at  the  moving  probe  with 
some  time  separation.  Because  turbulence  within  the  boundary  layer  is  nearly 
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homogeneous  in  x  and  convecting  downstream,  time  separation  can  be  used  to  estimate 
correlation  measurements  in  the  streamwise  direction  (Taylor’s  hypothesis). 

Because  of  the  physical  limitations  preventing  the  mapping  probe  from  reaching 
all  of  the  grid  points,  correlations  at  these  points  had  to  be  measured  using  additional 
experiments  or  by  interpolating  between  known  values.  Making  use  of  symmetry  in  z, 
some  missing  values  were  estimated  using  data  from  experiments  with  the  reference 
probe  at  other  locations.  For  locations  that  could  not  be  reached,  either  linear  or  para¬ 
bolic  interpolation,  depending  on  the  knowledge  of  the  correlation  variation  behavior 
in  space  and  time,  was  used  between  nearby  data  points  on  the  grid  to  estimate  the 
values. 

Case  2:  Outer  Laver  Structure  and  Large  Eddy  Break  Up  Devices 

The  boundary  layer  investigated  was  nominally  a  zero  pressure  gradient,  flat 
plate  boundary  layer  with  a  free-stream  velocity  of  10  m/s  and  Reynolds  number  based 
on  momentum  thickness  ranging  from  2000  to  8000.  A  flat  plate  manipulator  was 
placed  1.2  m  from  the  leading  edge  at  a  small  angle  of  attack  of  approximately  2 
degrees.  This  location,  x^,  corresponds  to  a  Reynolds  number  based  on  momentum 
thickness  Ree  of  about  2500.  The  LEBU  was  made  out  of  5  cm  wide,  0.127  mm  thick 
blued  spring  steel  strip.  It  was  suspended  1.9  cm  above  the  plate  and  was  carefully 
tensioned  to  avoid  any  vibration.  The  regular  boundary  layer  thickness,  Sq,  at  that  loca¬ 
tion  was  nominally  2.5  cm.  In  other  words,  the  LEBU  configuration  corresponded  to  a 
chord  length  of  2  Sg  and  a  height  of  0.75  Sg.  A  schematic  diagram  of  the  boundary 
layer  test  plate  and  of  the  LEBU  configuration  is  given  in  Figure  4.  For  more  details, 
the  reader  is  referred  to  Trigui  (1991). 

The  main  part  of  the  experiment  consisted  in  acquiring  six  simultaneously  sam¬ 
pled  hot-wire  signals:  three  velocity  components  at  a  fixed  reference  point  ("detection 
probe")  and  three  velocity  components  at  a  mapping  probe  traversed  over  a  coarse 
three-dimentional  sampling  grid.  In  order  to  study  the  evolution  of  the  large-scale  vort¬ 
ical  structures  as  they  pass  over  the  manipulator,  the  fixed  probe  was  positioned  0.635 
cm  (0.25  5g)  upstream  of  the  leading  edge  of  the  manipulator  at  the  same  height.  A 
mapping  was  performed  on  a  3x15x31  mesh  in  the  streamwise,  normal  and  span  wise 
directions,  respectively.  The  first  cross- stream  (y-z)  plane  corresponded  to  a  zero 
streamwise  separation  between  the  detection  and  the  mapping  probes.  The  second  and 
third  cross-stream  planes  coresponded  to  separation  of  1.27  cm  and  2.54  cm  (0.5  and 
1 .0  5g),  respectively. 
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Figure  4.  Schematic  Diagram  of  Boundary  Layer  Test  Plate  and  LEBU 
Configuration 
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Figure  5.  Schematic  Diagram  of  Data  Acquisition  Grid  for  Outer  Layer  Mappings 
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In  each  of  the  cross-stream  planes,  the  mapping  span  a  total  distance  ranging 
from  0.508  cm  to  3.81  cm  away  from  the  surface  of  the  plate  in  the  normal  direction, 
and  3.81  cm  on  either  side  of  the  test  plate  centerline  in  the  spanwise  direction.  These 
distances  correspond  to  distances  ranging  from  0.25o  to  I.SSq  in  the  normal  direction, 
and  from  -1.55o  to  I.SSq  in  the  spanwise  direction.  The  spatial  resolution  in  either  the 
normal  or  spanwise  direction  was  therefore  equal  to  0.254  cm  (0.1 5o).  As  mentioned 
above,  the  spatial  resolution  in  the  streamwise  direction  was  much  coarser,  namely 

l. 27  cm  (0.55o).  The  data  acquisition  grid  is  schematically  shown  in  Figure  5.  At 
each  grid  point,  long  time  records  of  the  six  velocity  components  were  acquired  at  a 
frequency  of  7,000  Hz,  corresponding  to  2  viscous  unit  time  units  between  samples. 
A  total  of  50,(XX)  samples  was  acquired  at  each  spatial  point.  The  digital  acquisition 
and  processing  procedures  were  very  similar  to  those  described  for  in  the  previous  sec¬ 
tion  for  the  Case  1  experiments. 

In  addition  to  all  the  usual  statistical  quantities,  such  as  means,  and  first  and 
second  order  moments,  all  nine  components  of  the  space-time  correlation  tensor  were 
computed,  yielding  /?„p  (Aa:,  Ay,  Az;  At)  in  both  the  regular  case  and  the  manipulated 
case.  From  this,  inference  about  the  underlying  spatio-temporal  structure  of  the  flows 
can  be  made,  and  other  quantities  such  as  convection  velocities,  integral  scales  in 
space  and  time  can  be  computed  in  the  usual  fashion.  Furthermore,  a  more  physical 
interpretation  of  these  extensive  correlation  measurements  was  made  by  using  the  sto¬ 
chastic  estimation  technique. 

m.  RESULTS 

Calibration  and  Use  of  Triple  Hot-Wire  Probes 

A  considerable  effort  was  devoted  to  establishing  a  convenient  and  reliable  cali¬ 
bration  for  arbitrary  geometry  triple  hot-wire  probes,  and  for  extracting  the  three  velo¬ 
city  components  from  anemometer  voltages  in  an  accurate  and  fast  manner.  The  con¬ 
ventional  models  to  represent  the  cooling  velocity  on  inclined  sensors  (Jorgensen, 
1971;  and  others)  were  extended  to  circumvent  some  of  the  known  limitations.  In 
addition,  a  novel  technique  was  designed  to  perform  the  inversion  of  the  non-linear 
calibration  equations,  regardless  of  the  form  of  such  calibration  equations.  The  pro¬ 
cedure  relies  on  a  clever  decomposition  of  the  three-dimensional  problem  into  a  series 
of  two-diemnsional  one,  which  are  then  solved  by  a  table  look-up  and  interpolation 
scheme.  This  new  technique  is  computationally  efficient  and  has  errors  which  are 
much  smaller  than  those  of  the  calibration  and  modelling  aspects  of  the  problem.  The 
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results  of  the  overall  procedure  were  presented  at  the  Symposium  on  Turbulence  at  the 
University  of  Rolla,  Missouri  in  September  1990  and  at  the  1990  APS  Meeting.  A 
paper  has  subsequently  been  submitted  and  accepted  for  publication  in  Experiments  in 
Fluids  (1991).  A  preprint  of  this  paper  is  included  in  Appendix  A,  which  describes  all 
the  details  of  the  technique. 

Case  1:  Wall  Layer  Structure  as  a  Function  of  Reynolds  Number 

The  bulk  of  this  aspect  of  the  project  is  described  in  details  in  the  Masters 
Thesis  of  Thomas  J.  Gieseke  (1990).  Unfortunately,  the  extensive  data  base  collected 
for  this  experiment  was  later  found  to  be  contaminated  by  cross-talk  among  the  hot 
wire  channels  and  has  not  yet  led  to  any  good  results.  The  cross-talk  problem  has 
been  identified  as  coming  from  the  use  of  the  sub-miniature  triple  hot-wire  probes, 
where  one  of  the  prong  is  common  to  ail  three  wires  to  minimize  the  probe  size.  We 
were  well  aware  of  the  possibility  of  cross-talk  at  the  time  of  the  experiments  and  had 
made  significant  changes  to  our  setup  and  anemometer  operations  to  eliminate  the 
problem.  Unfortunately,  some  cross-talk  remained  in  the  data  and  was  only  later 
identified,  thus  invalidating  the  data  as  it  stands.  Attempts  are  currently  being  made  to 
design  some  post-processing  schemes  to  eliminate  or  minimize  the  problem  a  pos¬ 
teriori.  We  have  also  ordered  some  new  subminiature  probes  with  no  common  prong 
(at  the  expense  of  a  slightly  larger  size).  These  probes  are  ordered  for  a  different  pro¬ 
ject,  but  we  expect  to  redo  the  boundary  layer  exjjeriment  in  the  coming  summer  to 
gather  a  new  data  set  (at  no  cost  to  the  current  project!).  We  believe  that  this  data  set 
will  yield  very  valuable  information,  particularly  in  light  of  the  progress  on  the  other 
aspects  of  the  project  (pseudo-dynamics  and  multipoint  conditional  averages). 

Case  2;  Outer  Laver  Structure  and  Large  Eddv  Break  Up  Devices 

A  very  extensive  study  of  the  outer  structure  of  the  turbulent  boundary  layer 
and  its  modification  by  Large  Eddy  break  Up  Devices  (LEBU’s)  was  carried  out  in 
conjunction  with  the  originally  proposed  tasks.  The  results  of  this  study  are  described 
in  details  in  the  PhD  Dissertation  of  N.  Trigui  (1991)  and  a  summary  of  the  some  of 
the  results  were  presented  at  the  1990  APS  Meeting  and  the  29th  Aerospace  Sciences 
Meeting  in  Reno,  Nevada  in  January  1991.  A  copy  of  the  AIAA  Paper  91-0519  is 
included  in  Appendix  B. 
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V.  CONCLUSIONS 

Based  on  our  composite  experience  from  the  various  studies,  it  has  become 
clear  that  the  stochastic  estimation  technique  is  an  extremely  valuable  tool  to  study  the 
structure  and  dynamics  of  turbulent  wall  bounded  flows.  The  "pseudo-dynamics" 
reconstruction,  coupled  with  the  use  of  multipoint  conditional  averages  provides  realis¬ 
tic  description  of  the  real  flow  structure,  otherwise  unavailable  and  masked  by  conven¬ 
tional  conditional  sampling  (single  condition  at  one  point).  The  richness  of  the  tur¬ 
bulence  and  the  dynamic  interaction  of  the  various  structures  is  completely  obliterated 
by  standard  technique  and  lead  an  erroneous  (or  at  least  very  incomplete)  picture  of 
the  flow.  The  results  obtained  thus  far  in  the  outer  region  of  the  boundary  layer  point 
to  a  succession  of  single-sided  eddies  arranged  loosely  in  a  staggered  pattern  in  space 
and  time,  very  similar  to  the  conceptual  picture  postulated  by  Stretch  and  Kim  based 
on  DNS  of  turbulent  channel  flow.  This  type  of  structural  organization  is  very 
different  from  the  conventional  symmetric  structure  picture  usually  discussed.  Conse¬ 
quently,  the  mechanisms  for  the  evolution  of  the  structures  may  be  entirely  different 
from  those  postulated  earlier  (for  example,  lift-off  away  from  the  wall  by  mutual 
interaction,  etc..). 

Furthermore,  we  have  recently  proposed  a  new  idea  (cf.  AFOSR  Workshop  on 
Turbulence  Structure  and  Control,  1991)  which  potentially  can  extract  real  dynamics 
from  the  stochastic  estimation  formalism,  alike  the  POD  and  dynamical  system 
approach  of  Aubry  et  al.  (1990).  Some  of  these  ideas  are  briefly  described  in  Appen¬ 
dix  D.  In  essence,  it  involves  using  as  "modes"  individual  sets  of  stochastic  estimation 
coefficients  from  a  multipoint  stochastic  estimate  and  plugging  those  "modes"  into  the 
Navier  Stokes  equations  to  derive  a  truncated  system  of  dynamical  equations.  Unlike 
POD  modes,  these  stochastic  estimation  "modes"  have  semi-compact  support  and  are 
not  unlike  an  optimal  wavelet  basis.  While  this  approach  is  very  speculative  at  this 
stage,  it  has  the  potential  for  major  breakthrough  in  our  understanding  of  wall  bounded 
flows  and  approximate  numerical  prediction  for  control  purposes. 
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NOMENCLATURE 


d 

Diameter  of  hot-wire 

Anemometer  voltage  output  for  the  t'*  hot-wire 

/ 

Length  of  hot-wire 

«,v  .w 

Orthogonal  components  of  velocity 

Dr\ 

Velocity  ratio  v/u  =  ranCv) 

Dr2 

Velocity  ratio  =  tan{^) 

Variable  ratio  Q^g/U^g 

Q 

Estimated  cooling  velocity  from  a  response  model 

Qi 

Cooling  velocity  associated  with  the  hot-wire 

Qmag 

Cooling  velocity  magnitude  iQi+Qi+Qj)'^ 

Qref 

Estimate  for  the  cooling  velocity  associated  with  the 
anemometer  output  at  a  known  probe  orientation 

Qr\ 

Cooling  velocity  ratio  j2i/G3 

Qr2 

Cooling  velocity  ratio  QAQi+Qh'^ 

RCd 

Reynolds  number  based  on  d  and  U^g 

U 

Velocity  vector 

Ub 

Velocity  component  normal  to  a  wire  and 
perpendicular  to  the  plane  of  the  wires  supports 
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U^g  Magnitude  of  the  velocity  vector 

Un  Velocity  component  normal  to  a  wire  and  in 

the  plane  of  the  wires  supports 
Ut  Velocity  component  tangent  to  a  wire 

e  Measure  of  error 

0,  Sweep  angle  for  the  i'*  hot-wire 

to,-  Angular  rotation  for  the  i'*  hot-wire 

(j),\|/  Angles  in  spherical  coordinates  describing  the  flow 

direction  with  respect  to  the  probe 

INTRODUCTION 

Hot-wire  anemometry  has  long  been  accepted  as  one  of  the  most  accurate 
and  affordable  methods  to  study  velocity  fields  in  turbulent  and  unsteady  Hows. 
Modem  experimental  methods  in  conjunction  with  powerful  computers  for  use 
in  data  reduction  have  allowed  for  the  development  of  probes  designed  to 
measure  the  three  components  of  velocity  concurrently  by  using  an  array  of 
three  hot-wires.  Each  wire  in  the  array  responds  to  the  cooling  effects  of  the 
flow  by  some  function  (2,  =  /,(«,v.w,0,  ,a>,)  where  0,  and  co,  describe  the  orien¬ 
tation  of  the  wire  with  respect  to  a  reference  coordinate  system  and  the  sub¬ 
script  i  denotes  a  specific  wire  in  the  array.  The  inversion  of  the  calibration 
equations  is  required  to  obtain  u  =  gi(Gi.i22-i23.6i.02-03-“i-®2.W3).  v  =^20.  etc. 
This  inversion  is  usually  complicated  and  tedious  because  the  cooling  equations 
are  coupled  and  nonlinear,  depend  on  probe  geometry,  and  contain  unknown 
constants.  Furthermore,  they  usually  yield  multiple  solutions  for  some  range  of 
cooling  velocities,  thereby  limiting  the  effective  flow  angles  which  can  be 
determined  with  a  given  probe. 

Few  ad-hoc  algorithms  have  been  developed  (Chang,  Adrian  and  Jones, 
1984;  Lekakis,  1988;  Vukoslavcevic,  Balint  and  Wallace,  1987;  Choi,  1989; 


Lekakis,  Adrian  and  Jones,  1989;  Russ  and  Simon,  1990;  and  others)  to  extract 
the  velocity  components  from  triple  hot-wire  measurements  and  most  rely  on 
certain  simplifying  assumptions.  The  most  crucial  assumption  is  the  choice  of 
the  model  to  describe  how  the  heat  transfer  from  each  wire  depends  on  the 
direction  and  magnitude  of  the  impinging  velocity  vector.  The  most  common 
model  takes  several  forms,  but  essentially  is  a  weighted  sum  of  the  velocities 
normal  and  tangential  to  the  wire.  The  probe  geometry  and  some  model  con¬ 
stants  must  be  experimentally  determined,  but  the  form  is  generally  accepted 
because  of  its  physical  basis  and  reasonable  agreement  with  experiment  (Mor¬ 
gan,  1975). 

The  resulting  set  of  equations  must  be  solved  point  by  point  for  the  velo¬ 
city  components  using  special  inversion  techniques.  Iterative  techniques  have 
been  used  successfully,  but  multiple  solutions  often  exist  and  the  solution  time 
is  lengthy  (Choi,  1989;  and  others).  Instead  of  computing  velocity  vectors  from 
hot-wire  measurements,  researchers  at  Stanford  University  (Yavuzkurt  et  at, 
1977)  developed  an  analog  circuit  to  solve  the  set  of  equations  in  real  time. 
Once  again,  the  governing  model  equations  had  to  be  simple  and  the  use  of  the 
circuits  was  limited  to  a  specific  arrangement  of  wires.  Lekakis,  Adrian,  and 
Jones  (1989)  developed  a  method  to  analytically  invert  the  three  cooling  veloci¬ 
ties  making  use  of  variable  transformations  very  similar  to  the  methods  used  in 
the  present  analysis.  The  inversion  procedure  requires  that  the  fomi  of  the  wire 
response  model,  or  cooling  law,  be  the  Jorgensen’s  form.  As  will  be  discussed 
later,  other  models  may  be  better  representations  of  the  cooling  law  when  prong 
interference  and  wire  wakes  have  significant  influences. 

Presented  here  is  a  set  of  procedures  designed  to  be  generally  applicable 
to  any  triple-wire  probe.  Wire  response  models  based  on  a  generalization  of 
the  physical  concepts  and  experimental  analysis  are  developed  for  each  indivi¬ 
dual  wire  within  the  array.  The  form  and  complexity  of  the.se  models  can  be 
arbitrarily  chosen  to  account  for  large  flow  angles  and  flow  interference 
between  wires.  A  method  to  efficiently  perform  the  numerical  inversion  of  the 
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calibration  equations  to  extract  the  velocity  vector  is  also  developed  which  can 
be  used  regardless  of  the  complexity  of  the  response  models  or  the  probe 
geometry. 

WIRE  RESPONSE  MODELING 

In  the  use  of  hot-wires,  the  wire  response  is  interpreted  using  the  concept 
of  effective  cooling  velocity,  Q.  The  cooling  velocity  associated  with  U  is 
defined  as  the  magnitude  of  a  velocity  vector  having  a  direction  normal  to  the 
wire  within  the  plane  defined  by  the  wire  supports  that  would  yield  the  same 
heat  transfer  from  the  wire  as  U .  It  is  assumed  that  the  anemometer  output 
voltage,  e,  depends  on  some  general  angular  response  model  and  on  the  wire 
sensitivity  relating  it  to  the  cooling  velocity.  We  will  show  that  the  probe 
response  (Q)  to  velocity  magnitude  (U^g)  and  velocity  direction  (U I U,^g)  can 
be  decoupled  by  determining  /  (U^g )  and  /  (U /U^g )  separately. 

Given  a  global  coordinate  system  defined  by  x  in  the  mean  flow  direction, 
y  in  the  vertical  direction,  and  z  in  the  spanwise  direction,  the  orthogonal  velo¬ 
city  components  associated  with  this  coordinate  system  («,  v,  w)  can  be 
transformed  into  a  local  coordinate  system  defined  by  the  wire  and  its  supports. 
The  two  geometric  parameters  required  to  define  the  new  coordinates  are  0,  the 
angle  at  which  the  wire  is  swept  back  from  an  orientation  perpendicular  to  the 
X  axis,  and  o),  the  angle  the  wire  is  rotated  about  the  x  axis,  measured  clock¬ 
wise  from  the  y  axis,  with  a  zero  degree  rotation  corresponding  to  an  orienta¬ 
tion  where  the  longer  of  the  two  supports  is  directly  below  the  shorter  one. 
The  new  coordinate  system  is  defined  with  one  axis  parallel  to  the  wire,  one 
normal  to  the  wire  in  the  plane  of  the  prongs,  and  the  third  one  (binomial) 
completing  an  orthogonal  set.  The  terms  U,,  U„,  and  Ui,  are  the  velocity  com¬ 
ponents  in  these  directions  as  shown  in  Figure  1.  In  terms  of  the  velocity  vec¬ 
tor  («,v,w)  and  the  wire  orientation  0  and  to,  these  velocity  components  are: 


(1) 
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t/^=  V  sino)  -  w  cosco 
U,=  u  sin0  -  V  COSO)  cos0  -  w  sinco  cos0  (2) 

Un  =  u  COS0  +  V  COSO)  sin0  +  w  sinco  sin0  (3) 

Experiments  on  heat  transfer  from  yawed  cylinders  (Morgan,  1975)  have 
yielded  approximations  for  the  cooling  velcxity  as  a  function  of  the  velocity 
terms  above.  The  most  commonly  accepted  form  is  Jorgensen’s  equation  (Jor¬ 
gensen,  1971): 

Q^  =  KiUSK2U^+K^U^  (4) 

The  values  of  Ki,  K2,  and  can  be  approximated  based  on  a  physical  argu¬ 
ment.  Insensitivity  to  tangential  velocities  implies  A")  =  0.  Equal  sensitivity  to 
the  normal  and  binormal  components  implies  K2  =  Kj=  1.0.  These  values  are 
commonly  assumed  because  they  greatly  simplify  equation  solutions  (Chang  et 
al,  1988;  Choi,  1989;  and  others).  It  is  usually  recognized  that  this  form  suffers 
from  several  shortcomings,  particularly  to  represent  data  at  large  flow  angles. 
However  due  to  its  simplicity  and  general  acceptance,  it  will  be  used  as  a  start¬ 
ing  point  for  our  hot-wire  response  analysis. 

It  has  been  found  (Balint  et  al,  1989)  that  the  "constants"  A",,  A'2  and  A'3 
should  be  allowed  to  vary  with  the  flow  angle  to  better  represent  data  at  large 
flow  angles.  An  alteinate  way  to  capture  these  effects  is  to  add  additional 
terms  in  the  Jorgensen’s  model  involving  higher  powers  of  U„,  and  U,. 
However,  those  higher  powers  of  the  velocity  components  must  be  appropri¬ 
ately  rescaled  by  the  velocity  magnitude  U^g  to  allow  the  separation  of  sensi¬ 
tivities  to  velocity  magnitude  and  flow  angle.  The  expanded  model  can  be 
summarized  as  follows: 

=  ATi  K2  -t-  AT,  + 


t// 


t/- 


+  a:. 


mag 


u. 


+  K, 


u: 


mag 


u, 


mag 


-!-  ••• 
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By  considering  a  spherical  coordinate  system  rather  than  a  rectangular  one, 
it  can  be  seen  that  the  three  velocity  components,  u ,  v  and  w ,  can  be  expressed 
in  terms  of  the  velocity  magnitude  U^g  times  functions  of  two  flow  angles 
with  respect  of  the  probe  axis  ^  and  \(/,  as  shown  in  Figure  2.  This  signifies 
that  the  velocity  magnitude  can  be  factored  out  of  the  three  velocity  com¬ 
ponents.  Similarly,  the  velocity  magnitude  can  also  be  factored  out  of  the  U„ , 
Ub  and  U,  components  since  those  are  linearly  related  to  m,  v  and  w  (Equations 
1-3).  Since  U^g  can  be  factored  out  of  the  three  velocity  components.  Equation 
5  can  be  divided  to  by  U^g  to  yield  a  right  hand  side  which  does  not  depend 
on  the  velocity  magnitude  but  only  on  the  flow  angles  $  and  y; 


=  /(<t>.V)  (6) 

This  separation  of  the  velocity  magnitude  from  the  directional  sensitivity 
allows  the  dissociation  of  the  probe  calibration  procedure  into  an  angular  cali¬ 
bration  test  and  a  sensitivity  determination  with  respect  to  the  velocity  magni¬ 
tude.  It  is  important  to  note  that  the  separation  is  allowed  by  the  form  of  the 
wire  response  model.  The  Jorgensen  equation  is  only  appropriate  for  flow 
angles  and  magnitudes  in  the  range  for  which  the  response  model  was 
designed.  A  typical  correlation  (Morgan,  1975)  indicates  a  similar  model  is 
accurate  for  U,/^[uJ+Ub  <  2.75,  Ud  =  200,  and  Reynolds  number  based  on  wire 
diameter  between  1  and  12.  When  using  a  5\xm  diameter  hot-wire,  these  Rey¬ 
nolds  numbers  correspond  to  velocity  magnitudes  between  3.4m/s  and  43.2m/s. 

Since  the  angular  response  is  primarily  governed  by  the  geometry  of  the 
probe  array,  the  Jorgensen’s  model  is  used  as  a  base  from  which  more  detailed 
response  models  are  determined.  These  more  detailed  models  are  used  to  cap¬ 
ture  effects  of  wake  and  prong  interferences  as  well  as  model  variation  with 
flow  angle.  The  response  model  will  not  change  over  time  so  the  full  response 
of  a  probe  in  terms  of  velocity  magnitude  and  angles  needs  to  be  performed 


only  once  (or  rarely).  With  the  angular  response  established,  the  probe  can  be 
calibrated  against  velocity  magnitude  at  a  fixed  orientation  (corresponding  to  its 
intended  use).  This  procedure  leads  to  considerable  simplifications  for  the  use 
and  calibration  of  such  triple  wire  probes  and  is  justifiable  as  long  as  the 
response  models  are  accurate  in  the  range  of  flow  angles  and  magnitudes  where 
the  probe  is  used. 

MODEL  EVALUATION  AND  TESTING 

In  order  to  determine  the  best  combination  of  variables  and  their  powers 
for  the  response  model,  experiments  were  conducted  to  establish  the  relation¬ 
ship  between  the  anemometer  output  voltage  and  the  flow  direction  and  magni¬ 
tude.  The  low  turbulence  wind  tunnel  at  The  Ohio  State  University  in  the 
Department  of  Mechanical  Engineering  was  used  with  a  special  positioning 
mechanism  mounted  in  its  test  section.  This  mechanism  consisted  of  a  rod 
with  one  end  fixed  in  space  and  the  other  sliding  inside  a  ball  joint  attached  to 
a  precision,  computer-controlled  traversing  mechanism.  This  allowed  the  rod  to 
be  oriented  at  any  angular  position  with  respect  the  free-stream  velocity.  A 
probe  holder  was  attached  to  the  rod  to  secure  the  probe  parallel  to  the  rod. 
The  probes  were  connected  to  Constant  Temperature  Anemometers  operating 
with  overheat  ratio  ranging  from  1.3  to  1.6  depending  on  the  probe.  The  out¬ 
put  of  each  anemometer  was  fed  to  a  King’s  law  type  analog  linearizer.  The 
purpose  of  such  circuits  was  not  to  linearize  the  hot-wire  signals  per  se,  but  to 
maximize  the  sensitivity  of  the  digital  acquisition  system  at  the  higher  flow 
velocities.  This  procedure  degrades  the  signal  somewhat  in  terms  of  the  noise 
contamination,  but  the  minimum  signal  to  noise  ratio  for  the  circuitry  was 
found  to  be  around  500  which  is  acceptable  for  our  needs. 

Several  triple-wire  probes  having  a  variety  of  wire  arrangements  were  stu¬ 
died.  Probes  1  and  2  (manufactured  by  Auspex)  were  of  identical  design  hav¬ 
ing  a  common  central  prong,  with  the  three  wires  swept  at  approximately  45 


degrees,  and  equally  spaced  120  degrees  from  each  other.  Probe  3  (TSI,  Model 
1295AV)  had  three  wires  swept  45  degrees  but  mounted  on  separate  prongs. 
Two  of  the  wires  were  arranged  in  one  plane  as  a  regular  X-wire  array  and  the 
third  wire  was  mounted  behind  the  other  wires  in  a  plane  perpendicular  to  the 
X-wire  plane.  Figure  3  depicts  the  geometry  of  this  probe.  For  the  remainder 
of  this  paper,  only  the  results  of  probe  3  will  be  shown.  Equivalent  or  better 
results  were  obtained  with  probes  1  and  2  for  reasons  that  will  be  discussed 
later. 

Before  attempting  to  determine  how  the  wires  in  the  probe  array  respond 
to  variations  in  the  flow  direction,  a  sensitivity  calibration  was  conducted.  A 
standard  experiment  to  determine  the  anemometer  circuitry  response  to  cooling 
velocity  is  to  mount  the  wire  in  a  flow  of  known  velocity  with  the  wire 
oriented  in  a  specific  position.  While  varying  the  speed,  a  data  base  can  be 
established  for  voltage  output  versus  velocity  for  one  specific  wire  orientation. 
This  yields  a  functional  dependence  between  the  anemometer  output  voltage,  e, 
and  a  rough  estimate  for  the  cooling  velocity,  j2r«/.  using  a  response  model 
with  assumed  constants. 

Qref=fie)  (7) 

The  resulting  calibration  curves  (G„/  versus  e)  served  as  a  reference  to  com¬ 
pare  wire  response  during  an  experiment  with  wire  response  during  the  angular 
calibration  tests.  As  expected,  the  curves  are  approximately  linear  due  to  the 
preconditioning  through  the  analog  linearizer.  A  cubic  polynomial  is  fit  to  the 
data  using  a  least  square  procedure  for  a  continuous  approximation  of  (2„/  over 
the  range  of  an  anemometer  outputs.  This  adds  to  the  uncertainty  in  measuring 
the  cooling  velocity  but  this  portion  of  the  experimental  approach  is  very  accu¬ 
rate.  Typically,  errors  in  the  curve  fit  are  less  than  0.25%. 

The  angular  calibration  was  conducted  by  varying  the  probe  orientation 
through  a  range  of  angles  with  respect  to  the  flow  and  the  velocity  magnitude. 
At  each  calibration  point,  the  free-stream  velocity  and  the  probe  orientation 
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were  converted  to  (m.v.w)  in  the  frame  of  reference  of  the  probe.  The  resulting 
sets  of  (u.v.w)  and  output  voltages  were  recorded  for  later  processing.  Typi¬ 
cally,  the  velocities  were  varied  over  a  three-dimensional  grid  defined  by  a 
range  of  flow  angles  of  ±  20°  in  the  horizontal  direction,  ±35°  in  the  vertical 
direction,  and  a  range  of  free-stream  velocities  commensurate  with  the  intended 
use  of  the  probes.  This  is  near  the  cone  of  acceptance  for  the  probes  in  certain 
directions.  Each  of  the  three  independent  variables  was  divided  into  at  least 
five  divisions  yielding  a  minimum  total  of  125  data  calibration  points.  The 
uncertainty  in  probe  positioning,  ±  0.4°,  was  found  to  be  the  most  significant 
source  of  error  in  this  calibration. 

The  angular  calibration  was  used  to  define  the  model  which  best 
represents  the  data.  As  a  first  approximation,  a  comparison  was  made  between 
Q,tf  and  Jorgensen’s  form  of  Q  with  all  model  constants  assumed  to  be  equal 
to  reasonable  values  (/fi=0.1,  and  a:3=1.0,  for  example).  This 

yielded  a  plot  of  Q  versus  j2r«/  as  shown  in  Figure  4  for  wire  3  (the  wire  most 
strongly  affected  by  prong  interference  and  wire  wakes).  If  the  sensitivity  cali¬ 
bration  and  model  were  exact,  there  would  be  a  linear  relation  between 
and  Q.  As  can  be  seen  in  Figure  4,  some  function  Q  =F{Q^^f)  can  be  fit  to 
the  data  points.  Since  the  points  almost  fall  on  a  straight  line,  a  quadratic  poly¬ 
nomial  was  deemed  sufficient  to  approximate  the  data  and  was  determined  by  a 
least  square  procedure.  The  average  error  is  1.6%  with  a  maximum  error  of 
7.9%.  It  can  be  noted  that  there  remains  a  significant  amount  of  scatter  around 
the  fitted  line.  This  can  be  attributed  to  the  fact  that  the  assumed  model  is  not 
exact.  To  improve  the  model,  Q  is  replaced  by  the  response  model  to  obtain: 

i/,2  Kj  i/„2  =  F\Q,,f)  =  A+  B  Q,,f+C  (8) 

One  constant  in  F^  or  must  be  fixed  in  order  to  determine  the  others.  Fix¬ 
ing  the  U„  coefficient  to  1.0  and  letting  all  the  others  vary  is  consistent  with 
the  concept  of  cooling  velocity.  Reordering  Equation  8  yields, 

=  i/,2  -H  ATj  i//  -  B  Q„f  -  C  ^ 


(9) 
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all  the  constants  in  which  can  be  determined  using  a  least  squares  procedure. 
For  wire  3  the  equation  was  found  to  be 

Q^=  0.063  U,^  +  1.03  +  1.00  = 

-121.0fll.5  a./ -  1-52^*/  (10) 

The  curve  Q  =F(Q„f)  is  shown  in  Figure  5a.  It  can  be  noted  that  the  scatter 
of  the  fit  has  been  considerably  reduced  from  Figure  4  by  finding  the  "best" 
constants  in  the  model  and  the  function  F.  The  average  error  has  been  reduced 
to  0.96%  with  a  maximum  of  5.0%.  Figure  5b  represents  the  data  obtained  for 
wire  1  and  is  more  typical  of  the  collapse  using  this  procedure.  It  has  an  aver¬ 
age  error  of  0.66%  with  a  maximum  of  2.8%.  Wire  3  (Figure '5a)  exhibits  the 
worst  behavior,  as  compared  to  other  wires  in  the  array,  due  to  its  inherent 
asymmetry  with  respect  to  spanwise  velocity  and  the  upstream  interference  of 
wires  1  and  2.  For  the  other  types  of  probes  investigated.  Figure  5b  is  more 
representative  of  the  response  of  all  the  wires. 

The  geometric  constants  describing  the  probe  geometry,  0,  and  o),  ,  have 
been  assumed  to  be  their  nominal  design  values  in  the  analysis  above.  A 
further  refinement  of  the  method  is  obtained  by  calculating  the  error  of  the 
curve  fit  Q  =F(Qref)  for  a  range  of  geometries  around  the  design  values  and 
picking  the  geometry  which  yields  the  smallest  error.  Typically,  the 
"improved"  values  of  the  wire  angles  specifying  the  geometry  of  the  probe  are 
found  to  be  within  a  few  degrees  for  their  design  values.  In  actuality,  they  may 
not  be  descriptive  of  the  actual  geometric  constants  of  the  probe  but  may 
represent  physical  aberrations  such  aerodynamic  prong  interference. 

As  Figure  5b  illustrates,  there  remains  little  scatter  in  the  experimental 
data.  Further  improvement  can  sometimes  be  obtained  by  adding  additional 
terms  in  the  model  as  stated  in  Equation  5.  In  the  case  of  our  probe,  the  third 
wire  is  mounted  behind  the  X-array  and  is  more  sensitive  to  a  positive  w  com¬ 
ponent  than  a  negative  one  due  to  the  inherent  asymmetry  of  the  design  with 
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respect  to  the  spanwise  direction.  For  the  third  wire,  the  additional  terms  in  the 
model  are  helpful  in  further  reducing  the  error.  On  the  other  hand,  the  addi¬ 
tional  terms  do  not  contribute  to  any  significant  improvement  for  wires  1  and  2 
(X-array).  The  response  model  for  wire  3  can  be  improved  by  adding  a  higher 
order  tangential  term  to  compensate  for  the  asymmetry.  This  is  equivalent  to 
Equation  5,  with  and  K(,  set  to  0.0.  As  shown  in  Figure  6,  the  collapse  of 
the  data  is  slightly  better  than  in  Figures  4  and  5a  having  an  average  error  of 
0.65%  and  a  maximum  of  4.4%. 


It  has  been  suggested  that  if  the  Jorgensen’s  response  model  is  used  in  its 
basic  form,  the  constants  AT 2,  and  are  not  really  constant.  As  noted  ear¬ 
lier,  the  variation  of  these  values  with  flow  direction  can  be  accounted  for  by 
the  addition  of  higher  order  terms  to  the  response  model.  For  example,  con¬ 
sider  the  case  where  the  tangential  component  of  velocity  had  increasing 
influence  on  the  heat  transfer  as  the  flow  became  more  tangent  to  the  wire. 
The  Jorgensen  response  equation  could  be  written  with  a  variable  U,  coefficient 
as. 


1.0  +  K„ 


u, 


mag 


(11) 


This  equation  is  mathematically  equivalent  to  Equation  5  with  different 
coefficients.  If  the  values  of  K’s  are  assumed  to  vary  with  velocity  magnitude, 
the  procedure  to  separate  angular  dependence  from  velocity  magnitude  depen¬ 
dence  no  longer  works.  Fortunately,  the  dependence  of  the  values  of  K’s  on 
the  velocity  is  small  over  the  range  of  velocities  considered  as  is  evident  from 
Figures  5  and  6. 


VELOCITY  EXTRACTION 


With  an  acceptable  response  model  established,  the  problem  now  becomes 
one  of  extracting  the  three  velocity  components  from  the  three  wire  outputs. 
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When  in  use,  each  wire  will  output  a  voltage  e  from  which  can  be 

computed.  It  then  remains  to  determine  U  by  inverting  the  three  model  equa¬ 
tions.  This  has  been  successfully  accomplished  in  the  past  using  special  pur¬ 
pose  hardware,  special  analytical  approaches,  and  iterative  techniques,  but  to 
analyze  the  amount  of  data  required  for  high  order  statistical  moments  on  large 
experimental  grids  accurately,  more  efficient  solution  techniques  are  required. 
One  such  approach  which  has  been  successfully  implemented  for  single  sen¬ 
sors  is  the  table  lookup  method.  In  principle,  this  approach  could  be  general¬ 
ized  to  the  three-dimensional  problem  by  using  a  set  of  three-dimensional 
tables.  However,  the  memory  requirements  to  store  a  finely  resolved  set  of 
tables  are  prohibitive  and  coarser  tables  give  rise  to  significant  interpolation 
errors.  In  an  effort  to  maximize  the  table  resolution  and  reduce  the  error,  a 
more  sophisticated  method  is  used  which  transforms  the  three-dimensional 
problem  into  a  two-dimensional  one.  This  approach  is  possible  since  the  form 
of  the  model  was  selected  such  that  the  velocity  magnitude  could  be  scaled  out, 
yielding  a  two-dimensional  problem  which  depends  only  on  the  flow  direction. 
Similar  variable  transformations  were  proposed  by  Lekakis  et  al.  (1989). 

All  discussions  that  follow  rely  on  the  assumption  that  the  response  model 
is  accurate  over  the  range  of  velocity  magnitudes  and  directions  in  which  the 
probe  is  to  be  used.  The  validity  of  this  assumption  is  based  on  the  accuracy 
of  the  model  determination  discussed  in  the  preceding  section. 

Taking  the  Jorgensen’s  model  as  an  example,  the  instantaneous  velocity 
magnitude  can  be  factored  out  to  yield: 

+  f(niUJUrnagf  (12) 

iQllUnu^gf  =  +  K22{U,IU^,f  +  K2,(UJU^^f  (13) 

iQ^'Umagf  =  +  K,2(U,IU^gf  +  (14) 

It  can  be  shown  by  geometrical  transformation  equations  that  the  right  hand 
side  of  each  equation  is  only  a  function  of  two  velocity  direction  variables 
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only.  By  taking  the  ratio  between  Equations  12  and  14,  a  new  variable  was 
defined  that  is  independent  of  velocity  magnitude, 

Q/fi  =  Gi/Q3  (15) 

Similarly,  two  additional  variables  can  be  defined  by: 

G/f2  =  G2/(G?-H23)’''  (16) 

Gm.,  =VG?-H2|+e3  (17) 

The  velocity  components  u,  v  and  w  can  be  transformed  in  a  similar  way  to 
yield  two  directional  variables  that  are  closely  related  to  the  flow  direction: 

D/jj  =  v/u  =  '  (18) 

D/i2  =  w/{u^+v^)^  =  tan(^)  (19) 

To  complete  the  set  of  equations,  the  velocity  magnitude  can  be  related  to  the 
magnitude  of  the  cooling  velocities  by: 

-  U^glQ^g  (20) 

This  equation  is  derived  by  summing  Equations  12,  13,  and  14  and  taking  the 
square  root.  Df^  is  defined  so  that  it  is  a  function  of  the  velocity  direction 
only,  as  were  the  other  variables. 

With  the  form  of  the  model  and  all  its  constants  determined  from  the 
angular  calibration  experiments,  a  data  base  can  be  generated  by  varying  the  U 
over  a  specified  range  commensurate  with  the  calibration  data  and  computing 
the  corresponding  Q  values.  This  data  base  is  determined  computationally  and 
made  as  fine  as  possible,  typically  using  176,400  data  points.  Applying  the 
appropriate  variable  transformations  described  above,  the  computational  data  is 
changed  to  the  three  variables  Dm,  D^,  and  as  a  function  of  Q^x  and  Qr2- 
These  three  variables  are  depicted  in  Figure  7a  through  7c  as  contours  plots  and 
gray  level  maps.  The  figures  can  be  thought  of  as  the  top  surface  of  a  three- 
dimensional  body.  Only  the  parts  of  the  surface  having  normals  with  com¬ 
ponents  out  of  the  page  are  considered.  The  ponions  of  the  body  concealed 
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behind  the  top  surface  represent  alternate  values  of  the  velocity  variables 
corresponding  to  given  values  of  Qj,  Q2,  and  Q3.  These  multiple  solutions 
occur  at  high  angles  of  the  flow  direction  where,  because  of  the  non-linear 
behavior  of  the  model  equations,  the  solution  manifolds  fold  back  on  them¬ 
selves  and  the  equations  yield  multiple  solutions. 

To  determine  this  cone  of  acceptance  of  the  probe,  an  experiment  can  be 
performed  where  a  probe  is  exposed  to  steadily  increasing  flow  angles,  and  the 
trajectory  in  the  variable  space  (G/fi.GRz)  traced.  The  path  begins  near  the 
center  of  the  shaded  space  and  progresses  toward  the  edge  of  the  region.  At 
some  flow  angle,  the  path  turns  around  and  proceeds  back  into  the  shaded 
region.  Where  the  path  reverses  is  the  boundary  of  the  region.  The  flow  angles 
that  correspond  the  boundary  defines  the  cone  of  acceptance.  It  can  also  be 
noted  that  the  boundary  of  that  domain  coincide  with  regions  where  the  isocon¬ 
tours  of  one  of  the  variables  become  very  close  together. 

The  location  of  the  velocity  variables  in  Qui-Qnz  parameter  space  can  be 
illustrated  by  plotting  points  with  coordinates  {Qrx^Qrz)  corresponding  to  all  the 
calibration  velocity  triplets.  This  type  of  plot  was  used  to  evaluate  the  probe 
performance  (see  the  section  on  errors).  Figure  8  shows  such  a  plot  with  the 
data  points  corresponding  to  the  uniform  distribution  of  m,  v,  and  w  from  the 
angular  calibration.  The  calibration  data  does  not  span  the  entire  cone  of 
acceptance  for  the  probe  which  is  represented  by  the  solid  line.  The  response 
models  were  applied  beyond  the  range  of  calibration  data  in  order  to  fully 
examine  the  equation  behavior.  For  this  reason,  the  outer  edges  of  the  figure 
represent  an  extrapolation  of  the  response  models  in  flow  regimes  where  they 
have  not  been  proven  to  be  accurate. 

Figure  9  shows  the  distribution  of  points  in  the  computed  data  base  which 
spans  the  entire  cone  of  acceptance  of  the  probe.  For  clarity,  only  a  small  sub¬ 
set  of  the  data  base  is  actually  plotted.  The  distribution  of  the  points  is  heavily 
skewed  to  one  side  as  a  result  of  the  probe  geometry.  The  location  of  the 
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velocity  point  (u,0,0)  is  found  at  (j2«i.G/?2)  ~  0.707).  The  most  limiting 

bounds  on  the  cone  of  acceptance  occur  when  the  flow  causes  large  or 
small  Qn  1.  The  effect  of  these  problems  is  to  limit  probe  use  to  flows  not  hav¬ 
ing  strong  V  and  w  components  and  to  reduce  probe  accuracy  for  certain  flow 
angles.  Probes  of  type  1  and  2  do  not  exhibit  this  strong  favoritism  to  one  side 
and  do  not  have  as  limiting  of  a  cone  of  acceptance  due  to  their  more  symmetr¬ 
ical  geometry. 

To  generate  the  tables  for  the  lookup  procedure,  the  following  procedure 
was  used.  The  data  base  was  sorted  by  (2/?i  and  and  a  uniform  grid  was 
overlaid  on  the  data  for  each  velocity  variable  as  shown  in  Fig¬ 

ure  9.  Since  the  computed  points  are  not  equally  spaced  in  this  parameter 
space,  an  interpolation  procedure  was  used  to  calculate  the  table  values  on  the 
regular  grid.  The  interpolated  values  were  obtained  by  a  weighted  average  of 
the  neighboring  points  closest  to  that  grid  location,  with  the  weights  defined  as 
a  function  of  the  distance  between  each  data  point  and  the  grid  point.  In  this 
fashion,  evenly  spaced  two-dimensional  tables  were  constructed  for  the  vari¬ 
ables  Dyji,  Dh2  and  D^3.  As  will  be  addressed  in  more  details  shortly,  very  lit¬ 
tle  error  is  introduced  by  this  interpolation  scheme  owing  to  the  density  of  the 
data  distribution. 

The  use  of  the  tables  is  strictly  limited  to  the  range  of  the  U  used  in  its 
generation  and  the  boundary  of  the  single  solution  domain.  Sets  of  cooling 
velocity  ratios  out  of  the  range  of  the  tables  yield  an  error  in  the  table  lookup 
which  must  be  accounted  for  on  an  od  hoc  basis  in  order  not  to  contaminate  the 
measurements  during  an  actual  experiment.  By  calibrating  the  probes  over  a 
sufficiently  large  range  of  flow  angles  for  the  intended  use,  such  errors  are  rare 
and  can  be  flagged  appropriately. 

To  use  the  tables,  the  following  extraction  procedure  is  used.  Measured 
voltages  are  converted  to  Qr*/  using  the  sensitivity  polynomial  calibration  of 
each  of  the  three  hot-wires.  These  are  then  converted  to  to  yield  the 
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estimated  cooling  velocities,  Q .  These  three  cooling  velocities  are  then 
transformed  to  the  variables  Qm,  Qr2^  and  Q^ag-  By  table  lookup  and  bilinear 
interpolation  between  table  values,  the  D^2  D/j3  velocity  components 
are  then  determined.  Transforming  these  velocity  variables  back  to  the  usual 


velocity  components  is  performed  by  a  simple 
Specifically,  these  equations  are  as  follows: 

coordinate  transformation. 

^/?3  Qmag 

ij  —  — ■  1  —  -  M  1 

(21) 

Vl  +  (1  +  ^R2  ^R\ 

V  =DRiU 

(22) 

w  =  £>/f2 

(23) 

Since  all  these  computational  steps  are  explicit  and  simple,  their  implementa¬ 
tion  is  computationally  very  efficient.  This  fact  is  the  chief  advantage  of  a 
table  look-up  procedure.  A  summary  flow  chart  of  the  overall  calibration  pro- 

cedure,  table  generation  and  velocity  extraction  is  given  in  Figure  10. 
ERRORS 


The  calibration  method  described  in  the  previous  paragraphs  model 
involves  assumptions,  interpolations,  and  curve  fits.  This  results  in  significant 
uncertainties  which  need  to  be  quantified.  Selecting  the  best  model  for  the  wire 
response  offers  the  most  flexibility  in  reducing  the  method  errors.  Figures  4,  5 
and  6  illustrate  how  these  errors  can  be  reduced  by  first  selecting  the  best  con¬ 
stants  for  the  basic  model,  and  then  by  adding  appropriate  terms  to  the  model 
when  it  is  warranted.  The  average  error  of  the  curve  fit,  in  terms  of  Q,  with  the 
assumed  constants  in  the  model  for  Figure  4  is  1.6%  of  full  scale  with  a  max¬ 
imum  error  of  7.86%.  As  mentioned  before,  these  errors  are  for  wire  3  (worst 
wire).  The  average  errors  for  the  other  wires  are  of  the  order  of  1%. 
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To  quantify  the  errors  in  terms  of  velocity,  the  following  procedure  was 
used.  For  each  data  point  in  the  angular  calibration  data  base,  the  velocity  vec¬ 
tor,  U ,  was  extracted  from  the  actual  anemometer  voltages  using  the  procedures 
of  the  previous  section.  Since  the  actual  velocities,  U^,  are  known  from  the 
probe  positioning,  a  measure  of  the  overall  error  can  be  defined: 


These  overall  errors  represent  the  errors  compounded  through  the  fitting  of  the 
model  and  its  inversion  through  the  table  lookup  procedure  and  interpolation. 
The  average  eiror  in  the  extracted  velocity  using  the  assumed  model  constants 
is  2.6%  with  a  maximum  of  12.1%.  By  finding  the  best  constants  for  the 
model,  the  average  error  for  wire  3  (Figure  5a)  can  be  reduced  to  0.96%  of  full 
scale  with  a  maximum  of  5.0%.  Similarly,  the  average  errors  for  the  other 
wires  are  reduced  to  values  of  the  order  of  0.7%  as  exemplified  by  the  excel¬ 
lent  collapse  of  the  data  in  Figure  5b  (wire  1).  Finally,  by  including  additional 
terms  in  the  model  for  wire  3  (Figure  6),  the  average  error  can  be  reduced  to 
0.7%  with  a  maximum  of  4.4%.  These  additional  terms  can  play  an  important 
role  when  the  range  of  flow  angles  is  large  for  the  geometry  of  the  probe.  The 
error  magnitudes  are  commensurate  with  the  error  usually  encountered  with  sin¬ 
gle  or  X-wire  probes. 

The  distribution  of  errors  is  shown  in  Figure  8  for  the  model  used  in  Fig¬ 
ure  5b.  Each  of  the  symbols  in  Figure  8  represents  the  velocity  error 
corresponding  to  an  individual  calibration  data  point.  It  can  be  seen  that  most 
of  the  errors  are  below  2%  except  in  a  few  regions  of  the  parameter  space. 

Finally,  a  potential  source  of  errors  stems  from  the  table  generation,  table 
look-up  and  interpolation  procedures.  To  investigate  this  type  of  errors,  a  com¬ 
putational  experiment  was  conducted.  Velocity  vectors  were  selected  over  the 
range  of  U  used  during  the  model  testing.  Using  the  response  model,  the 
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cooling  velocities  were  calculated  and  inverted  back  through  the  table  look-up 
procedure.  The  velocity  vector  determined  in  this  way,  U ,  was  compared  to  the 
initial  input,  Ua,  and  the  error  determined  as  before.  Through  the  numerical 
experiments,  it  was  found  that  tables  with  a  resolution  of  200x200  are  sufficient 
to  insure  an  accuracy  of  the  order  of  0.03%  for  the  average  error  with  a  max¬ 
imum  of  the  order  of  0.1%.  In  other  words,  the  numerical  errors  associated 
with  the  velocity  extraction  by  the  table  look-up  and  calibration  method  are 
negligible  with  respect  to  the  model  errors. 

APPLICATION  TO  MEASUREMENTS  IN 
TURBULENT  BOUNDARY  LAYERS 

These  calibration  and  extraction  procedures  were  tested  in  a  realistic  situa¬ 
tion  by  taking  data  in  a  turbulent  boundary  layer.  The  experiment  was  con¬ 
ducted  in  a  zero-pressure  gradient,  flat  plate  boundary  layer  at  a  Reynolds 
number  based  on  momentum  thickness  of  the  order  of  4000.  Long  records  of 
the  outputs  of  the  three  anemometers  were  digitally  acquired  at  a  sampling  fre¬ 
quency  of  2  viscous  time  units  =  2)  and  stored  to  disk  for  later  processing. 
This  was  repeated  for  a  range  of  heights  in  the  boundary  layer  in  order  to 
obtain  profiles  of  the  statistics  of  the  three  velocity  components. 

Figure  11a  depicts  the  mean  profiles  of  the  streamwise  (u),  normal  (v)  and 
spanwise  (w)  velocities  ploted  against  y/5,  where  6  is  the  99%  boundary  layer 
thickness.  As  expected,  the  streamwise  velocity  profile  is  in  good  agreement 
with  other  measurements  and  the  normal  and  spanwise  are  extremely  close  to 
zero.  It  should  be  noted  that  these  two  quantities  were  not  forced  to  zero,  but 
calculated  directly  from  the  data.  This  indicates  that  the  probe  calibration  per¬ 
forms  well,  including  in  regions  characterized  by  strong  spatial  velocity  gra¬ 
dients  and  turbulence  intensities.  Figure  11b  depicts  the  profiles  of  the  stream- 
wise  (u').  normal  (v')  and  spanwise  (w')  turbulence  intensities.  Again,  the  data 
is  in  good  agreement  with  results  previously  reported  in  the  literature.  Owing 
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to  the  size  of  the  probe  and  its  supports,  no  data  was  gathered  in  the  near  wall 
region.  The  profile  of  the  normalized  cross-correlation,  is  shown  in  Figure 
11c.  The  correlation  coefficient,  has  the  value  of  -0.41  throughout  the  log 
layer  which  is  as  expected. 

CONCLUSIONS 

In  conclusion,  an  experimental  approach  to  calibrate  triple  hot-wire  probes 
with  arbitrary  geometry  has  been  designed  and  tested.  This  method  provides  a 
way  to  optimize  the  constants  in  the  conventionally  accepted  model  and  to 
include  additional  terms  in  the  model  to  further  reduce  the  errors.  This  pro¬ 
cedure  decouples  the  angular  calibration  from  the  sensitivity  to  the  velocity 
magnitude,  in  a  manner  similar  to  that  used  by  Lekakis  et  al.  The  decoupling 
considerably  reduces  the  complexity  of  the  calibration  and  velocity  extraction 
procedures. 

By  scaling  out  the  velocity  magnitude,  the  three-dimensional  velocity 
extraction  problem  can  be  reduced  to  a  set  of  two-dimensional  problems,  which 
are  then  solved  by  a  table  look-up  procedure.  The  table  look-up  approach  has 
the  advantage  of  allowing  for  arbitrarily  complex  response  model  equations  to 
be  solved  provided  they  follow  the  rules  described.  This  was  found  to  be 
important  when  wire  and  prong  interference  changes  the  nature  of  the  flow 
within  a  probe  array.  This  method  is  computationally  efficient  and  does  not 
require  the  numerical  inversion  of  the  non-linear  model  equations. 
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Fifiurd  1  -  Definition  of  Global  and  Local  Coordinate  Systems 
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Figure  3  -  Geometry  of  Probe  3 
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Figure  4  -  Calibration  Curve  for  Wire  3  with  Assumed  Constants  in  Model 
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5  -  Calibration  Curve  with  Optimal  Constants  in  Model 
a)  Wire  3 


Q^efimls) 


b)  Wire  1 


Qimis) 


A-26 


Figure  6  -  Calibration  Curve  for  Wire  3  with  Additional  Terms  in  Model 


a  Function  of  Cooling  Velocity 
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Figure  8  -  Distribution  of  Calibration  Data  Points  as  a  Function  of  Cooling  Velocity 
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Figure  9  -  Table  Generation  Procedure 


ANGULAR  CALIBRATION  TABLE  GENERATION  VELOCITY  EXTRACTION 


Figure  10  -  Flow  Chart  of  Calibration,  Table  Generation  and  Velocity  Extraction  Pro¬ 
cedure 
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Abstract 

The  near  field  effect  of  a  single  Large  Eddy  Break  Up 
device  (LEBU)  on  the  oncoming  vortical  structure  in  a  tur¬ 
bulent  boundary  layer  was  investigated  in  great  details  by 
measuring  the  full  two-point  space-time  correlation  tensor. 
In  addition  to  conventional  turbulence  moments,  spatial  and 
temporal  integral  scales  were  determined.  The  stochastic  es¬ 
timation  technique  was  used  to  reconstruct  estimates  of  the 
full  three-dimensional  spatio-temporal  evolution  of  the  vorti¬ 
cal  structures  as  they  pass  over  the  manipulator  the  name 
"Large  Eddy  Break  Up"  device  is  indeed  justified  by  the 
results  obtained  here.  It  was  found  that  in  addition  to  the  ob 
vious  inhibition  of  the  normal  velocity  component  by  the 
LEBU,  a  number  of  other  interrelated  mechanisms  are  at 
play. 


Nomenclature 

Cj  Skin  friction  coefficient 

h  LEBU  height 

/  LEBU  chord  length 

RqP  Correlation  between  a  and  P 

Reg  Reynolds  number  based  on  momentum  thickness 

t  Time 

T  Integral  time  scale 

u  Streamwise  velocity 

v  Normal  velocity 

w  Spanwise  velocity 

X  Streamwise  coor^nate 

Xi,  Downstream  location  of  LEBU 

y  Wall  normal  coordinate 

z  Spanwise  coordinate 

Ftec-stream  velocity 
a  LEBU  angle  of  attack 

5  Boundary  layer  thickness 

52.  6  Momentum  thickness 

5^  Boundary  layer  thickness  at  LEBU  position 

Ar  Time  separation 

Ax  Streamwise  separation 

Az  Spanwise  separation 

^  Non-dimensional  distance  downstream  of  LEBU 

(x-Xi,)/8t 

Subscripts 

0  Refers  to  regular  case 
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Introduction 

The  effect  of  Large  Eddy  Break-Up  (LEBU)  devices  on 
turbulent  boundary  layers  has  been  studied  extensively  in  the 
past  decade.  These  studies  were  first  aimed  at  investigating 
the  potential  of  these  devices  for  viscous  drag  reduction. 
Among  these  studies,  one  can  cite  the  investigations  by 
Corke  et  al}^,  Hefner  et  al.  ^  Bertelrud  et  al.\  and  many 
others.  The  results  obtained  varied  considerably  from  one 
study  to  another  and  were  sometimes  conflicting.  While  there 
was  general  agreement  about  the  ability  of  these  devices  to 
produce  sizable  local  skin  friction  reduction,  there  was  con¬ 
siderable  controversy  to  the  extent,  if  any,  of  the  net  drag 
reduction  potential  of  these  devices.  This  prompted  a 
different  type  of  studies  aimed  at  optimizing  the  geometry  of 
the  LEBU’s,  as  the  parametric  investigation  of  Plesniak  and 
Nagib^.  These  investigations  indicated  the  marked  sensivity 
of  the  drag  reduction  results  to  the  geometry  of  the  LEBU’s 
and  identified  the  acceptable  parameters  to  achieve  drag 
reduction. 

However,  it  also  became  clear  that  a  mote  detruled 
understanding  of  the  effects  of  these  manipulators  on  the  tur¬ 
bulence  structure  was  requited  to  further  optimize  these  dev¬ 
ices  and  make  them  less  sensitive  to  minute  geometrical 
changes.  A  number  of  studies  dealt  with  documenting  the 
changes  in  the  various  turbulence  characteristics  downstream 
of  the  LEBU’s  and  identifying  some  possible  mechanisms  for 
these  changes.  Such  studies  have  been  conducted  by 
Guezennec  and  Nagib^,  Chang  and  Blackwelder^  Lemay  et 
al.*,  Trigui  and  Guezennec’,  and  others.  This  time  however, 
the  results  of  the  different  studies  were  in  good  agrement, 
namely  that  through  an  inhibition  of  the  entrainment  process, 
the  LEBU’s  reduces  the  scales  of  the  eddies  in  the  outer  and 
logarithmic  region,  which  results  in  a  drastic  decrease  of  the 
momentum  transport  and  therefore  a  decrease  in  the  local 
skin  friction.  This  conclusion  was  further  supported  by  the 
failure  of  the  LEBU’s  to  produce  any  sldn  friction  reduction 
in  fully  developed  turbulent  channel  flow  as  recently  reponed 
by  Prabhu  el  of.’®. 

Various  mechanisms  have  been  proposed  for  the 
modification  of  the  boundary  layn  by  the  L^U’s®,  such  as 
direct  inhibition  of  die  noniial  velocity  component  by  the 
manipulator,  and  the  interaction  of  the  shed  vorticity  in  the 
wake  of  the  LEBU  with  the  turbulent  boundary  layer  to  name 
a  few.  However,  these  mechanisms  are  primarily  inferred 
from  measurements  taken  downstream  of  the  LEBU’s  where 
it  is  hard  to  quantify  exaedy  the  relative  role  of  these  various 
processes.  While  the  direct  inhibition  of  the  normal  velocity 
component  by  the  LEBU’s  is  probably  the  most  dominant 
mechanism  in  the  immediate  vicinity  of  the  LEBU’s,  there  is 
no  detailed  documentation  of  the  near  field  of  the  manipula- 
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tor  and  the  alteration  of  the  vortical  structures.  A  more  fun¬ 
damental  understanding  of  these  mechanisms  is  essential  for 
future  design  of  manipulators  that  possess  stronger  and  longer 
lasting  effects  on  the  turbulence  structure.  Furthermore, 
unlike  the  initial  focus  of  the  earlier  work  on  drag  reduction, 
there  is  a  renewed  interest  in  these  devices  for  a  variety  of 
applications  ranging  from  local  skin  friction  reduction,  to 
alteration  of  the  transport  of  scalar  contaminants,  heat 
transfer  applications,  turbulence  suppression  for  optical  appli¬ 
cations,  boundary  layer  self-noise  and  drag  reduction  applica¬ 
tions  involving  a  combinations  of  techniques.  Experimental 
investigations  of  some  of  these  effects  have  been  conducted 
by  Trigui  and  Guezennec’’**,  Lindcmann’^  Sommer  and 
Petrie*^  Keith  et  al}*,  and  others. 

Up  to  date,  all  investigations  (except  one  flow  visuali¬ 
zation  study^^  to  the  authors  knowledge)  have  been  con¬ 
ducted  downstream  of  the  manipulators.  From  such  studies,  it 
is  very  difScult,  if  not  impossible,  to  unravel  the  detailed 
mechanisms  by  which  the  LEBU’s  actually  alter  the  tur¬ 
bulence  structure.  In  particular,  it  is  impossible  to  separate 
the  effects  of  the  manipulators  on  the  oncoming  large  scale 
structures  from  the  effects  of  the  shed  vorticity  on  what  is 
remaining  of  these  structures  downstream  of  the  manipula¬ 
tors. 
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This  paper  is  part  of  a  larger  study’*  in  which  the  near 
field  of  the  LEBU’s,  both  upstream  and  downstream,  was 
investigated  in  great  detail  to  elucidate  these  key  questions. 
This  paper  will  concentrate  on  the  effect  of  a  single  LEBU 
on  the  oncoming  large  scale  structures  in  the  boundary  layer 
before  it  reaches  the  trailing  edge  of  the  manipulator.  In 
addition  to  presenting  some  of  the  usual  one-  and  two-point 
statistical  turbulence  measurements,  a  detailed  documentation 
of  the  effects  of  the  LEBU  on  the  large-scale  voitical  struc¬ 
tures  will  be  presented  through  a  reconstruction  of 
conditionally-averaged  three-dimensional  velocity  field 
around  the  manipulators. 


Experimental  Procedure 

The  boundary  layer  investigated  was  nominally  a  zero 
pressure  gradient,  flat  plate  boundary  layer  with  a  free-stream 
velocity  of  10  mfs  and  Reynolds  number  based  on  momen¬ 
tum  thickness  ranging  from  2(XX)  to  80(X).  A  flat  plate  mani¬ 
pulator  was  placed  1.2  m  from  the  leading  edge  at  a  small 
an^e  of  attack  of  approximately  2  degrees.  This  location,  xi, 
corresponds  to  a  Reynolds  number  based  on  momentum 
thickness  Reg  of  about  2500.  The  LEBU  was  made  out  of  5 
cm  wide,  0.127  mm  thick  blued  spring  steel  strip.  It  was 
suspended  1.9  cm  above  the  plate  and  was  carefully  ten¬ 
sioned  to  avoid  any  vibration.  The  regular  boundary  layer 
thickness,  fig,  at  that  location  was  nominally  2.5  cm.  In  other 
words,  the  LEBU  configuration  corresponded  to  a  chord 
length  of  2  and  a  height  of  0.75  fig.  A  schematic  diagram 
of  the  boundary  layer  test  plate  and  of  the  LEBU 
configuration  is  given  in  IHgute  1.  For  more  details,  the 
reader  is  referred  to  Trigui’*. 

The  boundary  layer  characteristics  were  documented  by 
acquiring  a  set  of  velocity  profiles  on  the  plate  centerline  at 
various  downstream  stations.  The  location  of  the  various 


Fig.  1.  Schematic  Diagram  of  the  Boundary  Layer  Test 
Plate  and  LEBU  Configuration. 

measurment  stations  ranged  from  ^  =  0.5  to  53.0,  where  ^  is 
defined  as  the  non-dimensional  distance  downstream  of  the 
trailing  edge  of  the  LEBU.  A  total  of  49  points  was  collected 
in  each  profile.  The  points  were  equally  spaced  over  a  dis¬ 
tance  ranging  fix)m  100  wall  units  to  1.5  boundary  layer 
thickness  firom  the  wall.  At  each  point,  long  time  records  of 
the  three  velocity  components  was  acquired  with  a  triple 
hot-wire  probe  at  a  frequency  of  7,000  Hz.  A  total  of  50,000 
samples  was  acquired  in  each  record. 

The  main  part  of  the  experiment  consisted  in  acquiring 
six  simultaneously  sampled  hot-wire  signals:  three  velocity 
components  at  a  feed  reference  point  ("detection  probe")  and 
three  velocity  components  at  a  mapping  probe  traversed  over 
a  coarse  three-dimentional  sampling  grid  In  tnder  to  study 
the  evolution  of  the  large-scale  vortical  structures  as  they 
pass  over  the  manipulator,  the  feed  probe  was  positioned 
0.635  cm  (0.25  5g)  upstream  of  the  leading  edge  of  the  mani¬ 
pulator  at  the  same  height  A  mapping  was  performed  on  a 
3x15x31  mesh  in  the  streamwise,  normal  and  spanwise  direc¬ 
tions,  respectively.  The  first  cross-stream  (y-z)  plane 
corresponded  to  a  zero  streamwise  separation  between  the 
detection  and  the  mapping  probes,  '^e  second  and  third 
cross-stream  planes  coresponded  to  sqnration  of  1.27  cm  and 
2.54  cm  (0.5  and  1.0  5g).  respectively. 

In  each  of  the  cross-stream  planes,  the  mapping  span  a 
total  distance  ranging  from  0.508  cm  to  3.81  cm  away  from 
the  surface  of  the  plate  in  the  normal  direction,  and  3.81  cm 
on  either  side  of  the  test  plate  centerline  in  the  spanwise 
direction.  These  distances  correspond  to  distances  ranging 
from  0.25g  to  1.5Sg  in  the  normal  direction,  and  from  -1  .55o 
to  I.5Sg  in  the  spanwise  direction.  The  spatial  resolution  in 
either  the  normal  or  spanwise  direction  was  therefore  equal  to 
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Fig.  2.  Schematic  Diagram  of  the  Data  Acquisition  Mapping. 


0.2S4  cm  (0.15o).  As  mentioned  above,  the  spatial  resolution 
in  the  streamwise  direcdon  was  much  coarser,  namely  1.27 
cm  (OJSg).  The  data  acquisition  grid  is  schematically  shown 
in  Rgure  2.  At  each  grid  point,  long  time  records  of  the  six 
velocity  components  were  acquired  at  a  frequency  of  7,000 
Hz,  corresponding  to  2  viscous  unit  time  units  between  sam¬ 
ples.  A  total  of  50,000  samples  was  acquired  at  each  spatial 
point 

The  local  skin  friction  coefricient  along  the  plate  was 
evaluated  using  various  methods.  In  the  regular  boundary 
layer  case,  the  friction  velocity  was  determined  from  a  stan¬ 
dard  Qauser  plot  with  the  usual  constants  and  from  a  two- 
dimensional  momentum  balance  equation.  In  the  manipulated 
case,  since  it  has  been  shown  by  that  the  constants  in  the 
"law  of  the  wall"  do  not  hold  in  non-equilibrium  boundary 
layers,  only  the  two-dimensional  momentum  balance  was 
used  to  compute  the  local  skin  friction  coefficients  along  the 
plate.  The  practical  implementation  of  the  ntomentum  bal¬ 
ance  technique  was  similar  to  that  used  by  Corke  et  al}  and 
Plesniak  et  al?.  Although  such  technique  is  inherently  inac¬ 
curate,  it  is  the  only  tool  available  for  reducing  the  data  in 
manipulated,  non-equilibrium  turbulent  boundaiy  layers  short 
of  dfrect  measurement  by  local  skin  friction  balances’^.  In 
order  to  validate  the  use  of  the  two-dimensional  balance  tech¬ 
nique.  velocities  profiles  were  aquired  at  different  spanwise 
locations  on  either  side  of  the  plate  centerline  and  at  few 
downstream  locations.  It  was  found  that  the  boundary  layer 
thickness  did  not  vary  by  more  than  3%  over  the  span  of 


interest  Furthermore,  the  skin  friction  results  obtained  from 
this  technique  were  only  used  to  perform  a  relative  com¬ 
parison  between  the  regular  and  manipulated  cases,  hence 
minimizing  the  inherent  inaccuracies  of  the  technique. 

In  addition  to  all  the  usual  statistical  quantities,  such  as 
means,  and  first  and  second  order  moments,  all  nine  com¬ 
ponents  of  the  space-time  correlation  tensor  were  computed, 
yielding  /{qp  (Ax,  Ay,  Az;  Ar)  in  both  the  regular  case  and 
the  manipulated  case.  From  this,  inference  about  the  under¬ 
lying  spatio-temporal  structure  of  the  flows  can  be  made,  and 
other  quantities  such  as  convection  velocities,  integral  scales 
in  space  and  time  can  be  computed  in  the  usual  fashion. 
Furthermore,  a  more  physical  interpretation  of  these  extensive 
correlation  measurements  can  be  made  by  using  a  tool  called 
the  stochastic  estimation  technique.  This  technique  was  first 
introduced  by  Adrian  and  his  co-workers**’*®’®’.  Recently,  it 
has  been  us^  extensively  in  the  context  of  boundaiy  layers 
to  obtain  estimates  of  conditional  eddies  (in  tiie  sense  of  clas¬ 
sical  conditional  ensemble  averages)  by  extracting  informa¬ 
tion  from  the  (unconditional)  two-point  correlation  tensor. 
This  technique  is  very  powerful  and  convenient  when  dealing 
with  a  three-  or  four-dimensional  data  base,  particularly  to 
investigate  a  large  number  of  "conditions".  For  more  details 
about  the  derivation  and  the  use  of  the  technique,  the  reader 
is  referred  to  recent  applications  to  turbulent  boundary 
layers2'-22. 
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The  detailed  documentation  of  the  characteristics  of  the 
regular  and  manipulated  boundary  layers  can  be  found  else¬ 
where**.  The  results  are  similar  in  character  to  those 
reported  by  other  investigators  for  similar  cases.  As  an 
example,  the  downstream  evolution  of  the  skin  friction  ratio 
between  the  manipulated  and  regular  case  is  shown  in  Figure 
3.  A  maximum  of  about  20%  local  skin  friction  decrease  is 
obtained  and  the  effect  slowly  relaxes  up  to  ^  =  SO. 

A  comparison  of  the  downstream  evolution  of  the  vari¬ 
ous  integral  time  scales  is  shown  in  Figure  4.  The  left 
column  represents  the  integral  time  scale  associated  with  the 
streamwise  velocity  fluctuations,  while  those  associated  with 
the  normal  and  spanwise  velocity  fluctuations  are  shown  in 
the  center  and  right  columns,  respectively.  All  time  scales 
are  made  non-dimensional  with  the  local  mean  velocity  at 
that  height  and  the  momentum  thickness,  Sj.  In  effect,  these 
can  also  be  inteipreted  as  non-dimensional  length  scales  by 
the  use  of  Taylor’s  hypothesis.  The  normal  distance  from  the 
wall  is  also  made  non-dimensional  with  the  momentum  thick¬ 
ness.  A  decrease  of  about  20%  is  observed  for  the  integral 
scale  associated  with  the  streamwise  velocity,  T^,  and  this 
decrease  appears  to  persists  up  to  the  last  measurement  sta¬ 
tion  at  ^  =  S3.  On  the  other  hand,  the  integral  scales 
associated  with  the  normal  and  spanwise  velocity,  T,  and  T^, 
respectively,  are  very  significantly  reduced  in  the  upper  part 
of  the  boundary  layer  irruoediately  behind  the  LEBU.  This 
efrect  is  defimtely  less  pronounced  below  the  LEBU.  This 


-ao.oo  0.00  20.00  40.00  m.oo  10.00 


e 

Fig.  3.  Comparison  of  the  Downstream  Evolution  of  the 
Skin  Friction  Coefficient  between  Regular  and 
LEBU  case. 


indicates  that  different  mechanisms  are  at  play  above  and 
below  the  manipulator.  However,  both  T,  and  T,.  relax 
much  faster  than  T,  with  downstream  distance,  with  T, 
retuming  to  regular  values  by  ^  =  53.  These  relaxation 
effects  are  similar  to  those  reported  in  the  literature  for  the 
turbulence  intensity  profiles.  It  should  also  be  noted  that  all 
three  integral  scales  are  drastically  reduced  in  the  near  wake 
of  the  LEBU.  This  can  be  attributed  to  the  generation  of 
small  scales  in  the  near  wake  of  the  manipulator  which  dom¬ 
inate  the  boundary  layer  fluctuations. 

To  understand  what  is  happening  to  oncoming  large 
scale  structures  in  the  immediate  vicinity  of  the  LEBU,  all 


\  T.  u/^  T,  u/t, 


t'IS.S  <=IS.S  t.|S.S 


cjoB  fliHi  t.n  cn  i. 

L  u/f,  \  wtt  I.  wi, 


Fig.  4.  Comparison  of  the  Downstream  Evolution  of  the  Integral  Time  Scales  of  the  Three 
Velocity  Components  between  Regular  and  Manipulated  case:  a)  Streamwise  velocity, 
b)  Normal  Velocity  and  c)  Spanwise  Velocity. 
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nine  components  of  the  space-time  correlation  tensor  were 
measured  with  respect  to  a  fixed  point  located  just  upstream 
of  the  manipulator,  as  described  in  the  previous  section.  A 
comparison  between  the  regular  and  the  manipulated  case  is 
shown  in  Figures  5,  6  and  7  for  the  and  com¬ 

ponents,  respectively.  In  each  figure,  iso-correlation  contours 
arc  shown  in  the  r-z  plane.  The  time  axis  is  made  non- 
dimensional  with  and  5o  and  the  spanwise  axis  non- 

dimensionalized  by  Sq.  In  each  Figure,  the  left  column 
represent  the  regular  case  at  y/8o  =  0.6.  while  the  middle 
column  represents  the  LEBU  case  at  y/Sg  =  0.6  (just  below 
the  LEBU)  and  the  right  column  corresponds  to  the  LEBU 
case  at  y/Sg  =  0.9  (just  above  the  LEBU).  In  all  three 
figures,  the  top  line  represents  a  streamwise  separation 
Ax/Sg  =  0,  while  the  middle  line  corresponds  to  Ax/5g  =  O.S 
and  the  bottom  line  conesponds  to  Az/5g  =  I.O. 

As  expected.  Figure  5a  shows  that  Taylor’s  hypothesis 
holds  quite  nicely  for  the  streamwise  velocity  correlation 
in  the  regular  case.  The  convection  velocity  is  slightly  less 
than  as  expected  at  that  height  in  the  boundary  layer,  and 
the  structures  do  not  significantly  decay  in  space-time  over  a 
distance  of  a  boundary  layer  thickness.  A  similar  behavior  is 
observed  in  Figures  6a  and  7a  for  the  normal  and  spanwise 
velocity  correlations,  R„  and  R^,  respectively.  However, 
R„  does  decay  appreciably  more  than  the  other  two.  For  the 
manipulated  case  below  the  LEBU,  the  streamwise  velocity 
correlation  (Figure  5b)  indicates  a  similar  convection  velocity 
and  a  very  slight  strengthening  of  the  u  perturbation  as  it 


passes  under  the  LEBU.  This  can  possibly  be  attributed  to 
the  very  small  angle  of  attack  of  the  manipulator  contributing 
to  a  small  contraction  effect.  On  the  other  hand,  the  stream- 
wise  velocity  correlation  above  the  LEBU  (Figure  5c)  exhi¬ 
bits  a  convection  velocity  faster  than  probably  due  to  the 
acceleration  of  the  flow  above  the  leading  edge  of  the  mani¬ 
pulator.  Furthermore,  the  streamwise  velocity  perturbation 
decays  rather  strongly  in  space-time,  thereby  invalidating  the 
use  of  Taylor’s  hypothesis  even  over  shon  distance  in  the 
immediate  vicinity  of  the  LEBU.  The  difference  in  convec¬ 
tion  velocity  between  the  top  and  bottom  side  of  the  LEBU 
contributes  to  the  decorrelation  of  the  remnants  of  the  suuc- 
tures  over  the  manipulator,  ultimately  prohibiting  their  recon¬ 
nection  downstream  of  the  device.  This  difference  in  con¬ 
vection  velocity  is  clearly  linked  to  the  circulation  around  the 
device  and  partially  legitimizes  the  use  of  small  positive 
angles  of  attack  and  sufficient  long  chord  length^.  It  should 
also  be  noted  that  the  correlation  contours  are  slightly 
pinched  in  the  spanwise  direction  both  above  and  below  the 
LEBU,  Le.  the  manipulator  decreases  slightly  the  spanwise 
length  scale  of  the  structures^'^. 

On  the  other  hand,  the  normal  velocity  correlations, 

exhibits  similar  behavior  above  and  below  the  LEBU 
(Rgures  6b  and  6c).  The  structures  decay  very  rapidly, 
perhaps  more  so  on  the  top  side  Again,  this  points  to  the 
lack  of  validity  in  applying  Tay.or’s  hypothesis  even  over 
short  distances  in  the  inimediate  vicinity  of  the  manipulator. 
It  also  confirms  that  the  direct  inhibition  of  the  normal  velo- 
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Fig.  5.  Comparison  of  the  Space-Time  Correlation  Maps  of  the  Streamwise  Velocity  for: 
a)  Regular  Case  at  y/5  =  0.6,  b)  LEBU  Case  at  y/5  =  0.6  and  c)  LEBU  Case 
at  y  /5  =  0.9. 
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Fig.  7.  Comparison  of  the  Space-Time  Correlation  Maps  of  the  Spanwise  Velocity  for: 
a)  Regular  Case  at  y/8  =  0.6,  b)  LEBU  Case  at  y/5  =  0.6  and  c)  LEBU  Case 
at  y/5  =  0.9. 
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city  component  by  the  LEBU  is  felt  immediately  and  contri¬ 
butes  very  significantly  to  the  decay  of  the  vortical  structures 
over  distances  of  the  order  of  less  than  a  boundary  layer 
thickness  (less  than  half  a  chord  length  in  this  case). 

The  spanwise  velocity  correlations,  exhibit  a  very 
interesting  behavior  as  shown  in  Figures  7b  and  7c.  Not 
only  the  perturbation  does  not  remain  conelated  as  it  travels 
over  the  LEBU,  but  it  becomes  anti-correlated.  This 
behavior  can  probably  be  attributed  to  a  "splatting"  effect  or 
the  formation  of  a  new  vortical  structure  on  the  surface  of  the 
LEBU  due  to  the  newly  imposed  no-slip  boundary  condition. 
It  is  needless  to  say  that  the  concept  of  Taylor's  hypothesis 
loses  all  its  meaning  in  this  case. 

To  further  illustrate  the  structural  changes  to  the  vorti¬ 
cal  eddies  passing  over  the  LEBU.  the  stochastic  estimation 
technique  was  used  as  a  tool  to  reconstruct  estimates  of  the 


three-dimensional  conditional  velocity  field  associated  with  a 
typical  perturbation  impinging  on  the  leading  edge  of  the 
manipulator.  As  described  in  the  previous  section,  the  com¬ 
plete  space-time  correlation  tensor  is  used  as  a  basis  to  deter¬ 
mine  the  linear  stochastic  estimation  coefficients.  For  funher 
details  about  the  application  and  validity  of  the  stochastic 
estimation  technique  to  boundary  layers,  the  reader  is  referred 
to  the  work  of  Guezennec^*  arxl  Guezennec  and  Ttigui^. 
The  choice  of  the  "condition"  is  not  critical  to  the  conclu¬ 
sions  being  drawn  here.  As  an  example,  the  case  of  an 
irtrush  of  high-speed  potential  fluid  towards  the  wall  was 
chosen  as  the  relevant  "condition".  The  magnitude  of  u  was 
chosen  to  be  equal  to  one  r.m.s  value,  a„,  and  the  magnitude 
of  V  was  chosen  to  be  minus  one  r.m.s.  value,  a^.  The  mag¬ 
nitude  of  w  was  chosen  to  be  zero  (hence  yielding 
spanwise-symmetric  structures),  although  mote  likely  values 
be  randomly  ranging  within  plus  or  minus  one  r.m.s., 
and  would  yield  strongly  asymmetric  structures  as  pointed 


Regular  Case 

x/5  =  0.0  At  UJS  =  0.0 


LEBU  Case 

JC/8  =  0.0  At  Ujh  =  0.0 


Fig.  8.  Comparison  of  the  Space-Time  Evolution  of  Reconstructed  Velocity  Field  for: 
a)  Regular  Case  and  b)  LEBU  Case. 
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out  by  Choi  and  Guezennec^'*  (not  shown  here).  Figure  8 
shows  a  comparison  between  the  space-time  evolution  of 
such  a  conditional  strucmre  between  the  regular  case  (left 
column)  and  the  LEBU  case  (right  column).  The  top  line 
correspond  to  a  stream  wise  separation  Ax/8o  =  0.0,  i.e. 
upstream  of  the  manipulator  and  hence  both  cases  arc  identi¬ 
cal.  The  middle  line  represents  the  case  where  Ax/Sq  =  0.5 
and  A;  Ujbo  =  -0.5,  whereas  the  bottom  line  represents  the 
corresponding  case  for  Ax/Sg  =  1.0  and  A/  =  -1.0. 
Since  this  is  a  space-time  representation,  it  should  not  be 
interpreted  as  a  physical  spatial  structure,  but  rather  as  an 
illustration  of  the  concepts  discussed  in  the  previous  three 
Figures.  In  all  Figures,  the  cross-stream  plane  (y-z)  is 
shown  with  arrows  representing  the  v  and  w  velocity  com¬ 
ponents  in  that  plane.  The  gray  level  background  represents 
the  magnitude  of  the  streamwise  velocity  perturbation,  with 
white  corresponding  to  a  high  positive  perturbation,  nrid-gray 
to  a  zero  perturbation  and  black  to  a  high  negative  pemuba- 
tion.  The  scale  of  the  arrows  and  gray  levels  is  kept  fixed 
for  all  plots.  As  described  earlier,  Taylor’s  hypothesis  is 
quite  valid  in  the  regular  case,  with  the  structure  looking 
essentially  identical  for  all  three  space-time  locations  and  the 
amplitude  of  the  perturbation  slightly  decaying  for  u  and  w 
and  more  significantly  for  v.  On  the  contrary  in  the  LEBU 
case,  the  strength  of  the  vonical  motion  is  considerably 
weakened  by  the  the  firet  position  over  the  LEBU  and  the 
strength  of  the  u  perturbation  has  decayed  particularly  above 
the  LEBU  near  the  outer  edge  of  the  boundary  layer.  Again, 
it  should  be  re-emphasized  that  this  is  not  a  spatial 
representation  of  the  eddies  over  the  LEBU. 

Such  a  representation  is  given  in  Hgure  9,  which  dep¬ 
icts  a  comparison  between  the  regular  and  manipulated  case 
for  the  spatial  structure  with  zero  time  delay,  i.e.  a  true  spa¬ 
tial  structure  with  no  use  of  Taylor’s  hypothesis.  Again,  the 
"condition"  is  the  same  as  the  one  investigated  in  the  previ¬ 
ous  Figure.  The  regular  case  exhibits  a  pair  of  counter¬ 
rotating  eddies  which  are  inclined  over  the  wall  at  a  shallow 
angle.  By  the  last  station  (Ax/Sq  =  1.0),  the  eddies  are  not 
identifiable  since  they  have  reached  the  outer  edge  of  the 
bouitdaiy  layer.  In  the  LEBU  case,  we  can  observe  the  for¬ 
mation  of  a  secondary  pair  of  eddies  on  the  top  of  the  mani¬ 
pulator  at  Ax/Sg  =  0.5  due  to  the  "splatting"  effect  inferred 
earlia’.  The  remnant  of  the  original  eddy  passing  below  the 
manipulator  is  reasonably  unaHected  at  that  sution,  but 
represents  more  of  a  fossil  perturbation  not  associated  any 
more  with  a  vortical  structure.  By  the  last  station 
(Ax/Sg  =  1.0),  the  remaining  perturbation  is  considerably 
weaker  than  in  the  regular  case. 


Conclusions 

in  summary,  the  very  near  field  effects  of  a  single 
Large  Eddy  Break  Up  device  on  the  oncoming  vortical  struc¬ 
tures  in  a  turbulent  boundary  layer  was  investigated  with  very 
detailed  three-dimensional  measurements  of  all  velocity  com¬ 
ponents.  It  was  found  that  the  structure  of  the  boundary 
layer  is  affected  very  rapidly  and  in  many  ways  by  the  pres¬ 
ence  of  the  LEBU.  Since  all  three  components  of  the  velocity 
correlations  exhibit  a  different  behavior  as  the  vortical  large- 
scale  structures  pass  over  the  LEBU,  the  term  "Large  Eddy 
Break  Up"  device  can  be  better  understood.  The  mechanisms 


at  play  are  rather  complex,  ranging  from  the  "simple"  inhibi¬ 
tion  of  the  normal  velocity  component  of  the  eddies,  to  the 
formation  of  new  or  distorted  eddies  by  the  "splatting"  and 
the  no-slip  condition  on  the  manipulator  surface,  and  to  the 
shearing  by  different  convection  velocities  of  an  otherwise 
fossilized  streamwise  velocity  perturbation.  Furthermore,  all 
those  changes  happen  over  a  distance  of  the  order  of  a  boun¬ 
dary  layer  thickness  or  less.  While  this  does  not  elucidate  all 
aspects  of  the  role  of  the  LEBU  for  the  manipulation  of  tur¬ 
bulent  boundary  layers,  it  characterizes  in  a  very  detailed 
fashion  the  initial  effects  on  the  boundary  layer.  It  is  clear 
that  these  effects  arc  so  drastic  that  the  outer  structure  of  the 
boundary  layer  is  annihilated  for  practical  purposes  by  the 
time  the  turbulence  reaches  the  trailing  of  the  LEBU.  The 
additional  mechanisms  occurring  immediately  downstream  of 
the  manipulator  and  the  interaction  of  the  wake  of  the  LEBU 
with  the  boundary  layer  turbulence,  as  well  as  the  relaxation 
of  these  effects  with  downstream  distance  have  been  investi¬ 
gated  and  arc  described  elsewhere'*. 
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NEW  DIRECTIONS  FOR  RESEARCH  IN  TURBULENT  FLOWS 


•  Scientific  Needs 

Understand  routes  to  chaos,  transition 
— >  Understand  mechanisms  for  the  sustainment  of  turbulence 


•  Technological  Needs 

— >  Develop  better  transition/turbulence  models  for  flow  prediction 
-4  Develop  transition/tuibulence  control  schemes  to  achieve 
desired  characteristics 


•  A  few  facts  about  transition/turbulence 

Unsteady,  three-dimensional  flow  fields 
->  Complex  interaction  between  the  various  terms 
-4  Existence  and  importance  of  oi^anized/coherent  structures 
-4  Relative  lack  of  fonnalism  to  incorporate  heuristic  id^ 
into  models  or  dynamics 

-4  Lack  of  universality  among  flow  fields/researchers/models 
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IMPLICATIONS  FOR  RESEARCH  IN  TURBULENT  FLOWS 

•  Structures/flow  features  are  very  diverse.  Therefore  they  must  be  "captured"  in  the 
most  generic  way  to  avoid  bias. 

•  Structures  are  inherently  unsteady  and  three-dimensional.  Therefore,  the  methods  for 
their  "capture"  and  study  must  respect  those  attributes. 

•  We  must  incorporate  organized/coherent  structures  into  a  quantitative  formalism  to 
obtain  dynamics.  Therefore,  we  must  evolution  equations  for  these  stmctures. 

POSSIBLE  APPROACHES 

•  Stochastic  estimation  ->  Pseudo-dynamics,  p.d.f.  n^els,  L.E.S  subgm  models(?). 

•  Proper  Orthogonal  Decomposition  (POD)'— TD^ar^cal  systems  approach 
EXPERIMENTAL  NEEDS 

•  Measure  3-D  space(-time)  two-point  correlation  tensor  since  all  stractural  informa¬ 
tion  is  imbedded  in  it  Study  non-canonical  cases.  Use  either/both  approaches  to  ex¬ 
tract  specific  structures  and  study  them. 

THEORETICAL/NUMERICAL  NEEDS 

•  Use  experimental  data  to  study  the  following: 

-»  Three-dimensional  structure  (vclocityArorticity  fields) 

— >  Contributions  to  the  ovoall  flow  (kinetic  cno-gy,  Reynolds  stress,  produc¬ 
tion...) 

-4  Dynamics  (evolution  in  time,  regular/chaotic  behavior,  origin  of  chaotic 
behavior...) 

— >  Control  of  the  dynamics 

->  Sensitivity  of  results  (indirectly  of  correlation  tensor)  to  Reynolds  number, 
non-canonical  forcing... 


NECESSITY  TO  COMBINE 

EXPERIMENTAL/THEORETICAL/NUMERICALA^ISUALIZATION 
EFFORTS  TO  ENCOMPASS  THE  WHOLE  PROBLEM 


DRAWBACKS  OF  CONVENTIONAL  CONDITIONAL 
AVERAGING  TECHNIQUES 


•  The  conditional  average  decays  to  zero  away  from 
the  detection  point. 

•  The  choice  of  the  "condition"  is  difficult.  Results 
obtained  depend  directly  on  the  condition  selected. 

•  Strictly  kinematic  description  of  the  flow. 

•  Tedious  to  implement  for  many  different  "condi¬ 
tions"  for  a  3-D  data  base. 
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STOCHASTIC  ESTIMATION  OF  CONDITIONAL  AVERAGES 

-»  Least-square  estinudon  of  condkional  averages  (Adrian,  1976,  1979,  1980,  etc..)- 

(Conditional  average  of  u  at  x-H-  given  condition  at  x: 

0(x;r)  =  <  u(x;r)  |  u(x)  > 

Assume  a  linear  estimate  for  0(x:r): 

ii,(x;r)  =  Aijir)  Uj{x)  +  O  [«^(x)] 

Deteimine  Ajy(r)  by  least-square  procedure: 

*  <  [Afy(r)uy(x)  -  M,(x+r)J^  >  Min. 

Mmimize  the  error  with  respect  to  the  unknown  coefficients  A,y(r): 

^•(r)  ®  ®  *<>(x)(Aik  (r)ut  (x)  -  u,  (x+r)]  >  =  0 

which  yields: 

<  Uy(x)«j^(x)  >  A»(r)  =  <  My(x)u,  (x+r)  > 

The  deteimination  of  the  coefficient  A;y(r)  requires  the  use  of  the  2-point  2nd 
order  qMtial  coireladon  tensor. 

ESTIMATION  PROCEDURE 
-4  Measure  the  2-potnt  spatial  cotreladon  tensor. 

-»  Solve  linear  systems  of  equations  for  the  Afy(r)  coefficients  for  each  r. 

Select  condition  vectoifs)  at  x 

-4  Reconstruct  the  linear  estimate  of  the  conditional  velocity  field  around  x 
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GENERALIZATION  TO  INCLUDE  TIME  DEPENDENCE 

Conditional  average  of  u  at  x+r  and  /+t  given  condition  at  x  and  t : 

tl(x/  ;r,t)  =  <  u(x+r ;  /  -K)  |  u(x^ )  > 

Assume  a  linear  estimate  for  ti(x/  ;r,t): 

u,  (xj  ;r.t)  =  Aij  (r.t)  (x/ )  +  O  [u  \x^ )] 

Determine  Ajj{r,x)  by  least-square  procedure: 

e,-  =  <^Afy(r,T)u^(x/)  -  w,(x+r^-kc)j  ^  >  Min. 

Minimize  the  error  with  respect  to  the  unknown  coefficients  A;y(r,T): 

-  0 

dAij{r,x)  ~ 

which  yields: 

<  «y(x/)ut(x/)  >  Aa(r,T)  =  <  «y(x/)M,  (x+r/-K)  > 

The  determination  of  the  coefficient  A/jir.x)  requires  the  use  of  the  2-point  2nd 
order  space-time  correlation  tensor. 

ESTIMATION  PROCEDURE 
— ►  Measure  the  2-point  space-time  correlation  tensor. 

-¥  Solve  linear  systems  of  equations  for  the  A;y(r,T)  coefficients  for  each  r  and  each  t. 

— »  Select  condition  vectorfs)  at  x 

-♦  Reconstruct  the  linear  estimate  of  the  conditional  velocity  field  around  x  and  r. 
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GENERALIZATION  TO  IIICIIER  ORDER  ESTIMATES 

Assume  a  higher  order  expansion  for  the  estimate  of  u(x;r): 

(x;r)  =:  Af^(r)  Uj(x)  +  ujix)  u^ix)  +  Cfy*,(r)  «y(x)  «*(x)  «/(x)  +  •  •  • 

Determine  Afy(r),  BiJ^{r),  Qyufr), ...  by  least-square  procedure: 

e,- =  <  ^i«f(x;r)- «,-(x+r)j  >  Min. 

Minimize  the  error  with  respect  to  the  unknown  coefficients  Ajy(r).  Bij^{r),  Qyt/(r). ...: 

^iju{t‘) 

which  yields: 

<  tty(x)«*(x)  >  A/*(r)  +  <  My(x)«*(x)tt,(x)  >  Bij^(r)  + 

<  «y(x)(<t(x)u/(x)u„(x)  >  Cijuir)  -  <  Uj(x)uiix+xr)  > 

<  ay(x)«*(x)«,(x)  >  A,*(r)  +  <  «y(x)«*(x)«,(x)u^(x)  >  Bjj^ir)  + 

<  My(x)Ut(x)U/(x)U^(x)U.(x)  >  Cfyt,(r)=  <  Uy (x)^^ (x)U,- (X-hlT )  > 

etc... 

or  symbolically  represented  as: 

[AHCJ^IB] 

where  [A]  represents  the  matrix  of  moments  of  order  up  to  2"  at  the  reference  point  x. 
(C)  is  the  unknown  coefficients  in  the  stochastic  estimation  expression  (functions  of  r) 
and  represents  the  two-point,  n-th  order  spatial  correlations. 
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STOCHASTIC  ESTIMATION  OF  TWO-POINT  CONDITIONAL  AVERAGES 


-♦  Adritn,  1987. 


Conditional  average  of  u  at  r  given  condition  at  X|  and  X2: 


d(r;x,^2)  =  <  u(r)|  u(xi),u(x2)  > 


Assume  a  linear  estimate  for  u(r;X|4C2): 


«,(r;X|4C2) «y(xi)  +  B,y(r)  i<y(x2)  +  O  [a^l 


DetennineAjy(r)  and  by  least-square  procedure: 

€(  s  <[ifi(r) - ^  >  Min. 

Minimize  the  error  with  respect  to  the  unknown  coefficients  (r)  and 

dcj  dcj 

dAfy(r)  '  dBfjir) 

which  yields  symbolically: 

!A)IC]  =  IB1 


where  [A]  is  a  matrix  formed  by  the  second  order  one-  and  two-point  correlation  with 
point  X|  a^  X2,  [C]  represent  the  unknowns  coefficients  Ajj  and  Bjj,  and  [B]  rqKesent 
the  two^roim  qMtfial  correlation  between  pomts  r  and  X|  and  between  points  r  and  X2. 

•  If  the  two  reference  points  are  only  displaced  in  the  homogeneous  ditection(s) 
with  ttspect  to  each  ocher,  all  quantities  in  die  matrix  [A]  and  the  right-hand  side  [B]  can 
be  written  from  the  2-point  spatial  correlation  tensor /f;y(r;x). 

•  This  procedure  allows  to  defined  structures  of  a  qiecific  scale  (forced  by  the 
separation  be^een  X|  and  X2. 

•  Generalization  of  this  procedure  to  form  non-linear,  time-dependent  estimates  of 
multi-point  conditional  averages  can  similarly  be  obtained. 
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One-point  Condition 
t2  U/deltn  =  -.5 
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GOAL 


Design  a  tool  which  can  be  implemented  with  experimen¬ 
tal  data  to  examine  the  instantaneous  spatial  structure  of 
the  flow  and  its  evolution  in  time. 


PSEUDO-DYNAMICS  ^  REALISTICALLY  ANIMATED 

KINEMATICS 


N 
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MULTIPOINT  CONDITIONAL  AVERAGES 

•  Easy  to  implement  with  the  stochastic  estimation  foimalism  —  Nearly  impos¬ 
sible  to  implement  by  conventional  ensemble  averages. 

•  Control  over  timeAength  scale(s)  of  observation.  Study  the  interaction  of 
structures  rather  than  unrealisdc  structures  in  isolation. 


PSEUDO-DYNAMICS  REALISTICALLY  ANIMATED  KINEMATICS 

•  Choice  of  "conditions":  Measured  time  series  at  the  reference  points  to  get 
pseudo-dynamics,  i.e.  a  time  dependent  kinematical  representation  of  the  flow  con¬ 
sistent  with  a  measured  realization  of  the  flow  at  several  points. 


EXAMPLE  IN  THE  WAKE 
4  points  separated  in  time 


rONTnilB  LfUFlS 


High  stress  patch 


H’gh  stress  patch 


Plan  view 
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Estimate  Using  Wires  3  3c  6 

1000  data  points  \Mt* 


mV  4i*Vw 


o . G  Original 

* - *  Estimate 
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Wire  tfxariber 
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Estimate  Using  Wires  3  &  6 


1000  data  points 


Wire  Number 


V  Squared  Bar  U  Squared  Bar 
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